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RELAXING THE SUFFICIENCY CONDITIONS FOR ABSOLUTE STABILITY 


Vasile-Mikhai Popov 
(Power Institute of the Academy Science of the Rumanian People's Republic) 


The sufficiency conditions for absolute stability are investigated in the case of auto- 
matic control systems which contain a servomotor with a nonlinear speed character- 
istic. It is shown that it is possible to relax these conditions in certain cases, In par- 


ticular, necessary and sufficient conditions are obtained for three cases (except for 
intractable special cases). 


We consider an automatic control system described by the equations 


= >) + 41,..., 0), 
on n 


E=f(0), o= paye—t. 


where byg are constant parameters, 7 q are the coordinates of the controlled object, n, are constant parameters, 
€ is the coordinate of the controlling organ, p,, are controller parameters, £ (0) is a continuous function, hav- 
ing only a finite number of first-type discontinuities, and with the properties $(0) = 0 and o f(a) > O[class (A)). 


In [1] there is reproduced the proof of a theorem due to A. L. Lur‘e [2], determining the existence of criteria 
which, if satisfied, guarantee the absolute stability of a self-regulating system for any class (A) function f(a 
and, consequently, for any linear function f(o) = ha(0< h< @). 


Also in [1] (page 172) there was formulated the problem of inverting this theorem, i.e., the problem of 
ascertaining the validity of its converse, “that is, if a system is stable for an arbitrary function, f,, (0) = ho(0< h< o), 
then would it be correct to assert that the system would be stable for any class (A) function ?*. 


' This problem is solved in the case of three equations. The solution is based on a lemma, to be proved be- 


low, which allows the extension, in parameter space, of the domain of absolute stability, obtained beforehand by 
applying Lur'e's algorithm [2]. 


1. Auxiliary Lemma 


We consider an indirectly-controlled system described by the equations 


Zp =Apte +f (2) 4, 0), 
a= >) Bot 


with the usual assumptions (1, 2). 


| 


The following lemma is asserted: if the trivial solution to the system 


Zp = +f (3) (p= 1,2,..., 2), 


n (I) 


(where p is an arbitrary positive constant) possesses asymptotic absolute stability, and if this fact can be established 
by means of a Liapunov function of the form 


v= — O(2,, 2n)— | f (a)de, (1.1) 


(where @ is a positive definite quadratic form)and, moreover, if the following inequality holds: 


r+ (1.2) 


then the trivial solution of System (I) also possesses asymptotic absolute stability. 
We proved this lemma by constructing the Liapunov function for System (I). 


From the assumptions stated above, the total derivative of the Function (1.1), taking Equations (II) into ac- 
count, 


4(0) (3 Bote — rf (2)) — 2, (1.3) 
e=1 out 
is a positive definite function, 
We now prove that the function 


(1.4) 


_ which, by virtue of the assumptions already made, is negative definite, is indeed the Liapunov function for System 


(1). Indeed, the total derivative of this function, taking Equations (I) into account, has the form: 


= (6) (28,24 — (0) — 


(1.5) 


Since po f(a) >0 fora # 0, then (=) : turns out to be a positive definite function, q.e.d. 


It should be mentioned that the additional Condition (1.2) does not limit the generality of the argument, 
since this is a necessary condition for System (I) to be asymptotically absolutely stable ({2], page 75). 


If the parameters of System (II) satisfy the conditions previously established by Lur'e [1, 2] using Liapunov 
functions of the form (1.1), then the System (I) is asymptotically absolutely stable. It is proven below that, with 
a suitable choice of the arbitrary constant p, it becomes possible (except in certain special cases) to extend the 
region of absolute stability in parameter space. Let us carry out the analysis for the case h = 2. 


| 
| 
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2. Necessary Conditions for Asymptotic Absolute Stability 
We consider the system 


= ++ f (0) (p= 4,2), 
= + — rf (0). 


Replacing the nonlinear function £ (a) by the linear one 


1,(¢) = hs, (2,2) 
we obtain a linear system with the characteristic equation 


(2,1) 


+ + a, + + +5,) = 0, (2.3) 
where 


= —(Ay+ Ag), = Ary 
by = — (Ay + — (Bi + Be), = + Bade + Boda (2.4) 


For the roots of Equation (2.3) to have negative real parts, it is necessary and sufficient that the following 
inequalities hold: 


a,-+rh>0, 
hb, >0, (2.6) 
rh?b; + h(rdg + ab; — by) + > 0. (2.7) 
We assume that a;> 0 and a,> 0. Inequalities (2.5) - (2.7) will be satisfied for all positive h if and only 
if for 
r>0 (2.8) 
and 
b,>0 (2.9) 
the following alternative holds: 
by > Oand ray + a,b; — bg > —2 Vrajagd; , (2.10) 
or 
by = Oandra, > by. 


The last inequalities (2,10) determine the presence either of real negative, or complex, roots of the tri- 
nomial expression (2.7) in h, 


We rewrite Inequalities (2.10) in the following form: 


b, = Oandb, < ray. (2.11) 


Inequalities (2.8), (2.9) and (2,11) are necessary for the trivial solution of System (2.1) to possess asymptotic 
absolute stability. 


3 


! 
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3. New Fomulation of the Lur'e Conditions 


We consider the conditions for asymptotic absolute stability in the form [1, 2): 


r>0, (3.1) 
= — + + Bada 22 > 0 (3.3) 


(where € = 4 1) and we express Conditions (3.2) and (3.3) as functions of a, a2, by and b, (Cf. paragraph 2). Mak- 
ing use of Relationships (2.4), we get 
=2>0, (3.4) 
— [b, — Vra, Vb + ra, —a,b,) >0. (3.5) 


For Inequality (3.5) to be satisfied, it is necessary and sufficient that the roots of its left member, considered 
as a trinomial in Vb,, be real, and that the value of Vb, lie between these roots: 


b, >0, (3.6) 
e Vra, — < Vra, + Va,b;. (3.7) 


It is easily seen that the region corresponding to € = —1 is contained in that corresponding to € = +1, since 
Vb; > 0. Thus, the Lur'e conditions for ay > 0 and a, > 0 can be expressed in the following equivalent form: 


r>0, b, >0, b,>0, 
Vira, — < Vb, < + 


Comparing Inequalities (3.8) with the necessary Inequalities (2.8), (2.9) and (2.11), and excluding from con- 
sideration the boundary of the region of stability (for which by = 0), we are led to the conclusion that for 


a,b, > ra, (3.9) 
the Lur’e conditions are not only sufficient, but also necessary. If (3.9) is not satisfied, the Lur'e conditions no 
longer coincide with the conditions for necessity. * 


4. Extension of the Region of Absolute Stability 


We now apply our lemma to System (2.1). Repeating our previous reasoning (Cf. paragraphs 2 and 3) for 
the case of a system of the form (II) 


Zp = + (e=1,2), 


we obtain a characteristic equation of the form 


+ agh + h (rd? + + 04) = 0. (4.2) 


* It was already remarked in [1] that the boundary of the region obtained by means of Lur'e's algorithm coincides 
with the envelope of the family (2.7). 


| 
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Here, 


= + Bide + Badr + p + p 


where p is an arbitrary positive constant. 


Using the Lur'e conditions in the form (3.8), we obtain 


r,>0, b, >0, b, >0, 
Vra, — V < Vb, < + 


— ay, 


(4.5) 
(4.4) 


(4.5) 
(5.6 


By virtue of the lemma, if Conditions (4.5), (4.6) and (3.4) are met, the latter condition being replaceable 


by 
b,>0, 


since a,> 0, then System (2.1) is asymptotically absolutely stable. 


5. Decomposition of the Region D 


(4.7) 


As before (Cf. paragraph 3), we exluded fromconsideration the boundary of the region of stability, b, = 0, 


and we consider, instead of the regions (2,8), (2.9) and 2.11), the region 


r>0, b,>0, b,>0, 
Vb, <V ray + 


We decompose this region into three subregions; 


subregion Dy, where 
r>0, b,>0, b,>0, 
Vira, — < Vb, Vraz + Vayby; 


subregion D,, where 


r>0, b,>0, b,>0, 
Vb, Vayb;, — <0; 


and subregion Ds, where 


r>0, b, >0, b,>0, 
— > 0. 


It is easy to verify that the sum of the regions D,, D, and Ds equals 


D=D,UD,U Ds 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


(5.5) 


In region D, System (2.1) is asymptotically absolutely stable, by virtue of the Lur'e Conditions (3.8). Re- 


gions D, and Dy will be investigated below. 


6. Absolute Stability of System (2.1) in Region D, 


Indeed, in region D,, bj and bz (Cf. (4.3) and (4.4)] are positive for all positive p and, moreover, they are 


both monotone nondecreasing for increasing p. For p= 0, Inequality (4.6) does not hold in region D,. For increas- 
ing p, the leftmost member decreses monotonically to —a, the quantity Vb} remains monotone nondecreasing, and 
the rightmost quantity increases monotonically. Therefore, there is a positive value of P for which Inequality (4.6) 
holds. Since (4,7) also holds in region D,, all the conditions (4.5), (4.6) and (4.7) are satisfied in region D, and, 
consequently, System (2.1) is asymptotically absolutely stable in region D,. 


7. Absolute Stability of System (2.1) in Region Dy; 


In the region, bj is a monotonic increasing, and bj as montonic decreasing function of p. Therefore, it is 
necessary that the following inequality be satisfied: i 


P< — ath (7.1) 


Let us rid ourselves of the last constraint. 


If the following inequality is satisfied: 


(V ra, = V <b, <(V ra, +V ahi), (1.2) 


the Inequalities (4.6) and (7.1) are also satisfied, since from (7.2) it follows that bj> 0. 


If the inequalities r=20, bj > 0, b, > 0 and (7.2) are all satisfied, then the conditions established above (Cf. 


paragraph 4) are also satisfied, and therefore the trivial solution of System (2.1) possesses asymptotic absolute 
stability. 


Condition (7.2), having the form: 


M?<b,<N%, (7.3) 
is equivalent to the condition 
(b,)* — (M® + N*)b, + M*N?<0, (7.4) 
Making use of Relationships (4.3) and (4.4) and omitting the calculations, we obtain the inequality 
Ap* + Bp+C<0, (7.5) 
where 
A = (ra, — (7.6) 
4 B = ayb, (ra, — — (ra, + b,) (ray — by), (7.1) 
C= (ra, — a,b, (a,b, 2ra, — 2b,). (7.8) 
The trinomial on the left of Inequality (7.5) always has real roots in region Ds, since 
B* —4AC = 16ra,b, — + (74, (arb, — > 0. (7.9) 


In fact, in region (Cf. (5.4)] we have b, < ra, and b, > and consequently, 
ra,— b; > 0. (7,10) 


It remains to prove that in region Ds there exists a positive value of p for which Inequality (7.5) holds. 
We noticed that, for by = 0, the trinomial on the left of Inequality (7.5) has the multiple root 


(7.11) 


| 
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For increasing by, the roots of the trinomial become distinct, although remaining real [Cf. (7.9)]. Due to 
the continuous dependence of the roots on the coefficients of the trinomial, the roots will remain within the region 
p> 0 for sufficiently small by. If one of the roots of the trinomial in (7.5) becomes zero, then C = 0, The least 
value of by (Cf. (7.8)] for which this equality can bold will be 


= (Via — (1.12) 


but this is also the greatest value that by can assume in the region D,[Cf. (5.4)]. Consequently, in the region D, 
there is a positive value of p for which (7.5) is satisfied, and therefore System (2.1) is asymptotically absolutely 


stable. Since System (2.1) is asymptotically absolutely stable in regions Dy, D, and Ds, then by virtue of (5.5) we 
can formulate our conslusions, as detailed below. 


8. Necessary and Sufficient Conditions for Asymptotic Absolute Stability 


If it is assumed that ay > 0 and a, > 0, and if we exclude from consideration the boundary of the region of 
stability (by = 0), then the necessary and sufficient conditions for asymptotic absolute stability are: 


r>0, >0, b,>0, rag + Vayb;, (8.1) 


where by, by, ay and a, are the coefficients of the characteristic Equation (2.3) obtained from System (2.1) after 
linearization by means of the function f ,(0) = ho. 


Thus, it has been established that for three differential equations the lemma, proved in this article, allows 
one to obtain the region of asymptotic absolute convergence coinciding (with the exception of some special cases 


which require individual investigation) with the largest region of asymptotic absolute convergence which may be 
obtained in parameter space. 


The author desires to thank the investigators at the Institute Matematiki Akademii Nauk Rumynskoi Narodnoi 
Respubliki, working in the domain of differential equations, for their kind attention and valued comments. 
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OPTIMAL TRANSIENTS IN AN AUTOMATIC CONTROL SYSTEM HAVING 
A BOUNDED REGULATOR UNIT 


E. K. Krug and O. M. Minina 


(Moscow) 


The forms of the optimum transients in automatic control systems (ACS) 
are determined for those with objects of various response characteristics (includ- 
ing delays), assuming the control valve to be position bounded. It is shown that 
there are difficulties in using continuous-action regulators to produce optimal trans- 
ients, because the nonlinear converters to the regulators must have responses deter- 
mined by the magnitude and site of the perturbations, and by the initial values of 
the bounded co-ordinates. Stepwise operated regulators are proposed for use in 


obtaining optimal transients. 


In designing ACS one has to consider how to ensure the minimum duration of transients produced by per- 
turbations. The condition that the controlled quantity differ minimally from its set value is applied. 


_ To solve this problem one must analyze the behavior when the regulator is position bounded and various 
types of perturbation act, and produce transients of minimum duration. 


Object 


Regulator 


Fig. 1. Structural diagram of an ACS. 


Such transients will in future be termed optimal, 
the transient duration being taken as the time from when 
the perturbation is applied to when the controlled quantity 
shows precisely zero deviation. The latter situation can, 
as we know, occur in nonlinear systems. We must also 
determine the regulator response which provides optimal 
transients, and find ways of realizing it. 


Here we consider an ACS of block diagram as in 
Fig. 1. Here two types of perturbation act, one on the load 
(M) and the other on the input (X). The optimal transient 
responses will be found for stepwise perturbations to load 


and input for objects which are described by first-order differential equations with delays, and also for objects spec- 
ified by several first-order linkages. We suppose the regulator to have two extreme (bounded) positions. A fast 
regulator is assumed, i.e., one with time corstants incomparably smaller than those of the object. 


1. Derivation of the Optimal Transients in an Automatic Control System Having 


a Bounded Regulator Unit. 


In (1) we considered a linear ACS consisting of an object with a lag and a fast-acting regulator. It was 
shown that there was a limit to the transient duration even if the gain were infinite when there was a lag. The 
duration of the transient produced by a step perturbation to the load was 27, or r if to the input, r being the 
lag in the object (Fig. 2). The maximum deviation from the set value.in the output (y,,,,) when the perturba - 
tion is to the load can be determined if the response of the object and the magnitude of the perturbation are 
known. Thus for a first-order non-self -regulating object with a lag the maximum deviation (Fig. 2) can be found 


from 


| 
| 
| 
| 
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Ymax = AM (1) 
while for a first-order self-regulating object with a lag we have 
t 


= AMk(i—e (2) 


Ymax 


where AM= M~— Mp, Mg specifying the state before the perturbation, T being the time-constant of the object and 
k the gain in the object. 


X = const x 
M=const 
AX, 

t t 

Tt t 

+ t t 

Ln 


Fig. 2, wptimal transient response in an ACS of infinite gain. 


A regulator which is to provide the optimum transient response to such pertu:bations must, when the output 
deviates from the set value (y = 0), apply an infinite control signal and instantly set the controller in the posi- 
tion corresponding to the new equilibrium state, 


The fact that coordinates and their derivatives are always bounded in real systems is neglected in[{1]. Even 
if the regulator is very quick-acting and shows no coordinate bounding we cannot obtain the transients shown in 
Fig. 2 because the controller (throttle valve, slide valve, etc.) isposition-bounded, This corresponds to the ¢ co- 
ordinate (Fig. 1) being bounded. 


We shall derive the optimum transient response for the system shown schematically in Fig. 1 for a position- 
bounded controller. We shall assume the system in a steady state and that step perturbations to load and input are 
applied. It has been shown (2, 3] that the system will be of quickest response if the bounded coordinates have their 
maximum values during the transient. Hence in an ACS having an object with a first-order delay, the regulator 
must set the controller at the appropriate maximum displacement ft mg, When a perturbation acts at t= 0 to the 
input or att= 7 to the load, retain this setting for a definite time, and then set it at the position #,., correspond- 
ing to the steady state. 


In an ACS with such a regulator the maximum deviation y,,,,,. the transient duration t,, and the time the 
controller is held at the extreme position t* are related by the relationships given in Table 1 (jf, being the con- 
troller position before the perturbation). 


Table 1 is easily derived if we consider the transient in sections divided according to the controller position 
and assume the controller restored to the equilibrium position at a time r before the end of the transient. 


Fig. 3 shows the optimal transients for various step perturabations to the load for a nonself-balancing object 
with a first-order delay T/ r = 6 (Fig. 3a), or for perturbations to the input (Fig. 3b). 


The controller position bounding increases the duration of the optimal transient over that for an unbounded 
d _ system (Fig. 2), but the maximum deviation is unaltered. 


b 
M= const 


x, 


t t 


Fig. 3. Optimal transients in an ACS with a position-bounded 
controller, 


/ 


1 i 09 2 25X 


Fig. 4. Extreme position dwell time for controller vs perturbation: 
a) for a second-order object without self-regulation, b) for a second-order object with 
self-regulation , c) for a third-order object without self-regulation. 


a 
| 
Ms Xs 
M, X, 
| M; 43 
Xx 
4, 
Mg t t 
M, 
| M, x; 
T. t Tt 
My 
| M, 
} “2 2 
| t, t, ty 
| 
| 2 4 h 2 
| 
10 


The transients shown in Fig. 3 and Table 1 relate to alternate perturbations to a system initially at equili- 


brium (y = 0, dy/ dt = 0, # = Mg). 


TABLE 1 


Parameter 


Object with a first-order delay 
without self-regulation 


Object with a first-order delay 
with self-regulation 


y for 0<t<t 
y for 


y for 2n<t<t, 


ss 


yfor 0<t<t 


yfor t<t< 


Step perturbation to load 
y=0 


y= (t—*) 


Xx (t — 2t) 


AM 
Ymax~ T * 
2(Ymax— Ho) — 4M 


max — #0 — 


AM = 
Vmax — 4M 


c= 


=H + AM 


If perturbations are applied to load and input the relationships of t*, 


¢ 


y=0 
y=AMk(i—e 7) 


¥max= 4Mk 


4AM 
2t4+T — op x 


x F) 44] 


AM 
* 


44] 


ss = Mot AM 


Step perturbation to input 


y= AX 


max 
AX 


Ho 

Ax 


*+T in 


t* =Tin 


if both act together and are steps the extreme position dwell time can be found from 


AXT 


v= 


max — 4M 


If the load step occurs at a time r before the input step we should use 


tp and Ymax Will be different. E. g-, 


AXT + 
Umax — 4M 


From the above we may draw up the requirements for a regulator which will provide the optimum transient 
in a position-bounded controller system containing a first-order object with a delay. 


(4) 


Y= Ht may — | 
— AMe 
— Yo — 4M) 
¥max 
i 
= 
y= Ax 
) 
Ymax Ymax 
AXT 
|| = Ho 
11 


1, When the controlled quantity deviates from its set value the regulator must instantly set the controller 
in one of the limiting positions, if both perturbations act together. 


2. The controller must be held at the extreme position for a definite time t* , depending on the object 
parameters and the type of perturbation. 


3. The controller must then be set at u.., corresponding to the new steady state in the system. 
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Fig. 5. Transients with various types of controller reversals. 
Full lines: a) two, b) three reversals; dashed lines, one reversal 
(ty = t® = XT/ 


A regulator which fulfills these functions can only produce optimal transients in an ACS containing a first- 
order object with a delay, We will determine what transients we will get with such a controller in an ACS with 
a second-order object, and how they differ from the optimal ones. 


The optimal transients in an ACS containing an n-th order object and a position-bounded controller can be 
obtained by switching the controller n times from one limiting position to the other. Fig. 4 gives calculated 
curves for the relation of controller extreme position dwell time to input signal for various objects. It is clear 
that if one of the time-constants in the object is greater than the others the transient duration (t,+ t,+ ts+ ...) 
is mainly determined by the dwell time for the first position (t). 


Fig. 5 shows the optimal transients in systems with second- and third-order objects, using numerical ex- 
amples. Transients for the same ACS with one controller reversal are also shown, the controller dwell time being 
found from the expression for t* (Table 1), by replacing T by the largest time-constant, and 7 by the sum of the 
others (noting that t*< ty). 
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The examples of Fig. 5 enable us to derive the deviations from the optimal transients in the latter (one re - 
versal) instance, for high-order objects. The examples given in [4] indicate that the actual transients with one 
reversal are close to optimal for their high-order objects. Hence we may assume that the regulator problem for 
such objects is very similar to that for first-order objects. 


The regulator problem may be solved in two ways. 


1. By inserting nonlinear converters in a continuous-acting regulator [isodrome (proportional + floating), 
derivative, etc.]. This regulator must ensure that the controller is set in definite positions at set times and pro- 
vide a stable steady state. 


2. By producing a discontinuous-acting regulator with the same functions, 


2. Inserting Nonlinear Converters in a Continuous-Acting Regulator to Obtain 
Optimal Transients. 


As it is difficult to produce signals proportional to the second, third, etc., derivatives of the controlled co- 
ordinate in real apparatus, the simplest continuous-acting regulators have been widely used in practice (isodrome, 
proportional, first-derivative). Such systems are analyzed and approved from linear studies [5, 6}, the nonlineari- 
ties occurring in real systems being neglected. 


It has been shown [7] that controller position bounding increases the transient duration and deviation greatly 
with large perturbations. 


Fig. 6. Transients for an ACS with a first-order object and a delay: 

a) without self-balancing, perturbation to load, b) the same, perturbation to 
input, c) with self-balancing perturbation to load, d) the same, perturbation 
to input. 


Fig. 6. shows the optimal transients for an ACS with a first-order object with a delay (dashed) plus those 
in an ACS with an isodrome regulator of paiameters selected in accordance with Ziegler and Nichols" rules [5] 
(full lines). Fig. 7 shows how t, varies with the size of the perturbation for an ACS with a first-order object (non- 
self-balancing) and a delay.* It is clear that the transient duration is reduced by a large factor if the regu- 


. tp for an ACS with an isodrome regulator is considered as the time from when the perturbation is applied to 
when y = 0.01 X,. 
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lator provides optimal transients. 


We may wonder whether the transients may be brought near to optimal by using nonlinear converters in con- 


tinuous-acting regulators. 


The nonlinear relationships required for this purpose can be found by considering the behavior during optimal 
transients in the (dy/ dt, y) and (y, Jydt) planes. This is possible because the controller has definite fixed positions 


16 

72 

Perturbation 

10 to load 
8 

6 

4 

2 

4M 4X 
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Fig. 7. Relation of transient duration 

to magnitude of perturbation. Full lines: 
optimal ACS without coordinate bound- 
ings dashed: optimal ACS with coordi- 
nate bounding; broken lines: ACS with 
an isodrome regulator allowing for posi- 
tion-bounding in the controller. 
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Table 2 summarizes the equations describing the 
phase behavior in these planes for the above objects in re- 
lation to magnitude of perturabation. Figs. 8 and 9 give 
the phase trajectories for various perturbations. The ini- 
tial conditions for input perturbations are taken as 


y=AX, dy/dt=0, 


For load perturbations (Fig. 9) the initial conditions are 
y=0, dy/dt=AMIT, p=po. 


The relations between dy/ dt and y or between y 
and f ydt in Figs. 8 and 9 for which the controller exchanges 
an extreme position for the equilibrium one are termed 
the equations of theswitchover lines (see Table 2). These 
are easily found as the points on the trajectories from which 
the image point reaches the equilibrium position in a time 
t for a given perturbation. 


Figs. 8 and 9 and Table 2 show that these equations 
are not the same as those for the phase trajectories. The 
switchover equations differ according to the magnitude 
and nature of the perturbation (whether to load or input). 


For the controller to return from the extreme posi- 
tions at defined instants we must have nonlinear converters 
with characteristics corresponding to the switchover lines. 
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Fig, 8. Phase pictures for input perturbations. 
The dashed lines are for switchovers when X x X,4 (curve 1). 


The case X> — > X, corresponds to curve 2 in Fig. 8a and point 2 in 
in 8b, In Figs. Ba and b 0, ~ Ho 
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If, for instance, we use a first-derivative regulator with control law 


(y + /at), 


(5) 


and if the gain k is large enough (as soon as y= 0 the controller moves to the extreme position), then for a load 
perturabation AM > 0.5 (may ~ Mo) the regulator will hold the controller fully over for a time t = t* such as to 


give the optimal transient. 


TABLE 2 


Equations 


Parts of phase picture 


in the plane J vat; v 


in the plane y 


Perturbation to load 


Ph for T 
tajecto- 
es max 
» —AM< vit 
ver 
4(p max Ho) 
— max— #0) Try. — Ty?) 
Perturbation to input 
Phase for 0<t<t 
ydt = AXt + — 
=7 (Umax — 
for 4X> 
>F max 


dy AM trax + He 


[dy / dt}? 


“max 


= [dy / dt}, 


y= Ax 


t 


max 


= (Umax — Ho) 


When 4M < 0.5(H may ~ Ho), Or for an input perturbation, we cannot get the corresponding optimal transient 
with such a regulator, since the relationship between the proportional and derivative terms in (5) will not correspond 


to the switchover lines (see Table 2). 


Fig. 10 shows how y, the control signal y + rdy/ dt and @ vary during the transient when a step load perturba- 
tion occurs (fy = 0, AM= 0.75 max) O8 the condition that the regulator coefficients are related as in (5), Fig. 
10 shows that the control signal changes sign at t= r + t®*, the controller transferring from +#/-,,, t0~#max» and 
not to H.. = AM. This results in self-oscillation. It is thus clear that even if we could (by knowing the type and 
size of perturbation and initial controller position) reset the controller from the extreme position at a given mo- 
ment, we cannot thereby get it back to the steady state; this results in self-oscillation. 


If the switchover functions are generated by using nonlinear converters in the isodrome (proportional + floating) 
regulators the resulting system is structurally unstable. In fact, if the control signal, consisting of proportional and 
integral terms, is to be zero at some set time the integral term must have the reverse sign, 


These peculiarities, particularly the fact that the nonlinear converters are not single-valued, which is re- 
lated to the differing initial conditions, and the instability in the equilibrium state, do not only occur in systems 


with delays. 


Let us consider a system with a second-order object and a position-bounded controller, To obtain optimal 
transients here we must have two times (t, and t,, Fig. 6b) when the controller is reset to extreme positions. 
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Fig. 9. Phase pictures for load perturbations. 
The dashed lines are for switchovers when AM < 0.5 imax (curve 1) 
and when 4M > 0.54 may, (curve 2). 
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Fig. 10. Transients in an ACS with non~self-regulating first-order 
object with a delay. 


Fig. 11 shows the phase trajectories corresponding to the optimal transients in response to input perturbations 
for two initial states (fy = 0 and fy = 0.75 Hay) in dy/ dt and y coordinates. The switchover line (dashed in Fig. 
11) is not single-valued, and also depends on the perturbation and initial state. 


If the controller here switches over at given times but does not reach the equilibrium position at t= t+ t, 
we shall again have self-oscillation. 


To produced optimal transients in an ACS the continuous-acting regulators plus nonlinear converters must 
contain elements which effect the switchover functions as appropriate to the magnitude, type and site of the per- 
turbation and to the controller position. They must also contain elements which enable one to set the controller 
in a position corresponding to the perturbations applied. This rules out the regulator operating continuously dur- 
ing the transient. 


However, this does not imply that we cannot improve the response of an ACS subject to perturbations of a 
single species by using nonlinear links. 


3. ACS with a Discrete-Action Regulator 


Here we describe the operation ofa discrete -action regulator with which we can obtain optimal transients 
in an ACS containing a first-order object with a delay, and also consider examples of the transients resulting from 
various perturbations. 
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Fig. 11. Phase trajectories for an optimal ACS with a second-order object. 


The conclusions on the optimal transient in response to load perturbations are drawn on the assumption that 
the magnitudes of these are not measured directly; but the expressions of Section 1 ( (3) - (4), Table 1) contain 
the load perturbation. 


If the load perturbation can be measured we can design a regulator which can balance out the perturbation, 
without interference from the lag in the object. The transient duration will then be determined only by the regula- 
tor lag, and systems containing objects with very different parameters and quick-acting regulators will only show 
transients in response to input perturbations. 


In an ACS of object response described by linear equations we can evaluate the load changes (if the object 
patameters are known) from the way the controlled coordinate and its derivatives change. For first-order objects 
with delays the derivative of the deviation specifies the load change in the time preceding the derivative determi- 
nation (AM= Tdy/ dt). In such cases, by measuring the controlled coordinate and its derivatives we can design 
a discrete -action regulator having the functions stated inSection 1. 


In an ACS of this type t* and #,, can be found from suitable formulas using a computing unit. 


The computer must be supplied with y and dy/ dt at defined times. These times must not correspond to ex- 
treme positions in the controller. With objects having delays it is advantageous to measure y and dy/ dt at a time 
r after the controller has been adjusted from an extreme to a calculated position. 


The requirements as to the functions of a discrete -action regulator are thus as follows: if the ACS was at 
equilibrium (y = 0, dy/dt= 0, = pp) and the output deviates because of some perturbations, then the regulator, 
having measured some y, and dy,/ dt, must set the controller in an extreme position, hold it there for a time tf 
and then set it at the new position #,. At time 7 after setting the controller to fy the regulator must measure 
fresh values y, and dy,/ dt (again test the system), and from these calculate t,*; if t,* differs from zero 
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the controller is set in an extreme position for time t,* and then at pf, for a time r, etc. 


If t,* is zero the controller must pass directly to 2. At atime f¢ later y and dy/ dt are again determined, 
etc. 


If tz = O and py= py the ACS is at balance. In this state y and dy/ dt are measured at intervals of r. 


b 
t t 
d 
t 
t 


Fig. 12. Transients in an ACS with a discrete-action regulator. 


If any additic ial perturbation occurs between tests (perturbations applied subject to arbitrary initial condi- 
tions) the system wil react next time y and dy/ dt are tested. 


In this ACS the discrete-action regulator will sample y and dy/ dt at intervals of r+ t*. If the perturbations 
are small the intervals between tests will be close to r. 


An ACS with a first-order non-self-balancing object with a delay has relations which define t* and p,, of 
the form 


(y+ tdy/dt)T 


ss = Pins = 4 + Tdy/ dt (i=0, 1, 2, 3,...). (7) 


Figs. 12 and 13 show the transients in an ACS with the above discrete-action regulator for step perturbations 
of different sizes simultaneously applied (Fig. 12a), or applied at different instants (Fig. 12, b and c), for uniformly 
changing perturbations to the input (Fig. 12d), and for sinusoidal perturbations to the input (Fig. 13). 


Table 3 gives the conditions on which Fig. 12 was drawn up. 


Thus a regulator of this type can provide optimal transients to step perturbations to load and input and to 
uniformly changing input perturbations. 


We would expect such a regulator to give near-optimal transients with other slowly varying perturbations, 
since the interval between tests depends on the deviation in the output and on its derivative. 
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Since this regulator sets the controller in a position corresponding to the steady state for a time ¢ between 
samples while allowing for the perturbations applied, it provides stable operation. 


Fig. 13. Transients for sinusoidal input perturbations with > = 6, x m = 9-4(H may ~ Ho)- 
The full lines are for an ACS with a discrete- ~action reguiatee, the dashed lines for one 
with an isodrome regulator. 


We also note that using computers can substantially improve the response in systems in which the time taken 
to displace the effector is comparable with the time-constants of the controlled object. 


TABLE 3 
Fig. |Curve Perturbation characteristics 
i2,a | 12 M = MAndy X, for t<0, M = M,+ AM andX¥=X,+AX for #>0 
1 AM=0.5 (umax— Ho), AX = 0.25 (u max— He) 
2 4M =0.5 — Ho), AX =0.75 (4 max— He) 
12, 1,2,3 M = M. <0, M=M. f , f 
b | 1,2 for t<0 ot+-AM for 
1 AX =0.25 (Umax — Mo) AM = 0.25 (may — Ho) 
2 AX =0.25 (tmax — tHe)» AM = 0.5 (vmax— He) 
8 AX =0.25 (p max ~~ He)» AM = 0.75 — Ho) 
12, c | 1,2,3 X=X,for t<0, X=X,+AX for for for 
1 AM = 0.5 max — Ho)» AX = 0.25 (4 max— Ho) 
2 AM = 0.5(umax — AX =0.5 (& may— Ho) 
12, 4 | 2,2,3 M=M,—const, X= X, for 
1 Ax =0, k,=0,1 
max Ho 
2 AX = 0.07 
— Ho 
3 AX = 0.75 (max — He)» k, 


The relations defining t*, which allow for speed bounding in the effector, are easily derived for step per- 
turbations to load and input. 
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Relations giving t® for first-order self-regulating objects with delays, and for the objects of higher order con- 
sidered at the end of Section 1, are also easily derived. In the latter case the expressions for t* will contain not 
only the deviation and its first derivative, but also second, third, etc. derivatives. 


SUMMARY 


1. If optimal transients are produced in an ACS with a position-bounded controller the transient duration 
can be greatly reduced and the dynamic accuracy increased, as compared with ordinary systems. 


2. It is difficult to use continuous-action regulators with nonlinear converters to provide optimal transients 
in practice because the nonlinear links must vary with the type, magnitude and site of the perturbation and with 
the initial values of the bounded coordinates. Self-oscillation can also occur in such systems. 


3. The proposed discrete-action (sampling) regulator can provide optimal transients in response both to step 
perturbations to load and input, and to uniformly changing input perturbations. The formulas for the controller 
extreme position dwell time (6) and the equilibrium values of controller coordinates (7) are derived for this 
regulator. 
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EVALUATION OF THE ACCURACY OF INPUT REPRODUCTION FOR 
LINEAR SERVOS AND RECORDING SYSTEMS 


G. Vasil'ev 
(Moscow) 


The necessary and sufficient conditions for a linear system to reproduce a 
given class of inputs accurately are given. The maximum moduli of the inputs 
and their local rates of change are used to define the class of inputs. 


The paper extends earlier ideas by others [1-5]. 


1. General Concepts and Terminology 


1, Reproducing system, By this is meant a device of any operating principle which is intended to transform 
an input signal into an output, subject to a present relation between their time relationships. 


Such systems are subdivided into simple reproducing, differentiating, and integrating, depending on the type 
of correspondence between input and output. 


_ For a simple system the relation of output &,,,, (t) to input a; ,, (t) is described either by 
Gout Ka (t), (1) 
= which expresses an idealized followup process, or by 
Sour (0 = (2) 


which expresses an idealized recording process. 


A system which produces the desired correspondence absolutely precisely will be termed a precise reproducing 
system. Hence one which effects (1) will be a precise servo system of scale factor K, while one which effects (2) 
will be a precise recording system of scale factor K and time displacement r. 


Strictly speaking it is impossible to produce a precise system to operate on arbitrary inputs in practice. The 
reason is that there are various sources of reproduction error. 


Here we consider the reproducing properties of linear systems with lumped constant parameters. The repro- 
duction errors in such systems arise from their selective reactions to the inputs. 


2. Transfer and distortion functions, The transfer functions for systems of this class belong to the class of 
rational functions of the complex variable p. 


In general, the transfer function of a linear reproducing system K(p), can be represented as the quotient of 
two polynomials in ps 


1+ ap 


To this correspond a precise servo system of transfer function 


Ky(p)=K(0) (K@=lim K(p)) (4) 
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and a precise recording one of transfer function 


1) 
K.(p) = K (0) e"* (= tim (5) 


It is clear that the natures assumed for Kp) and K,(p) are such that they can be found exactly by calcula- 
tion or experiment. 


In fact K(0) appearing in (4) and (5) is numerically equal to the output produced by a real system from the 
input 
ain(t)=1(t) 00), 


and KOoy K(0) is numerically equal to the slope of the phase-frequency characteristic g(w) at w = 0. 


It is also advantageous to use some accessory functions (distortion pee to specify the reproduction pro- 
perties of the system, namely: 


Xo(p) = (6) 


K )— K, (p 
(1) 


By using the distortion function we can consider the real system as two independent systems acting in parallel, 
k()(0) 
one being a precise system of transfer function Kop) = K(0) (or A.(p)= K(0)e * _ ) and the other a distortion 


system of transfer function Agp) = K(0)x [or A,(p) = K(0)x 


3. The input. The law followed by the input is usually unknown. Certain aspects of the class of inputs are 
frequently known, however. 


The first is that all inputs in the class satisfy the Dirichlet conditions in the interval 0 <t < o, plus the 
condition 


for 0>t>—oo. 


If these conditions are fulfilled it is possible to effect a Laplace transform for each of the inputs: 
Ain (p) = \ ain(t) 


Most frequently, however, the properties of the input are expressed via more concrete and practically con- 
venient formulas. 


Thus in the usual frequency-spectrum interpretation® the properties are expressed in terms of the maximum 
frequency of the spectral components (in that class) and the limiting (for that class) modulus of the inputs: 


M= t)}. 
max sup |4iq(4)| (8) 


* This interpretation assumes that one can effect a Fourier transform 


Ajn (i) = (t) at 


for any of the probable inputs in the class. 
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Some [6-8] have used the basic functions of random process theory (correlation and spectral density) to spe- 
cify the characteristics of the class of inputs. 


In the above cases the suitability of the system for reproducing a given class of inputs is analyzed via the 
frequency characteristics, which illustrate the behavior of the transfer function along the imaginary axis in the P 
plane. 


Other methods are possible. Thus in[2, 3] a complex-spectral interpretation of the input is used;* in ad- 
dition to M above, the coordinates and dimensions of the region where the poles in the Laplace representation 
of the probable inputs in the p plane lie, are used. 


Here we use M and Q to express the properties of a given class of inputs. 
Parameter Q is given by 
a (t) 


and is henceforth termed the limiting local rate of change index (for that class); it applies to the time interval 
0st=o.°** 


Q=max sup sup 
0<k<o 


(9) 


Values of Q are given below for certain elementary types of input. 


Type of input sin wt wt 


Corresponding value of Q 0 a o o 


As an accessory means of representing a given class of inputs we use the configuration and coordinates of 
the region where the poles in the transforms lie. This region is henceforth denoted by R, , and the corresponding 
class by the By class, specified by M and Q. 


The system suitable for reproducing a Pn class can be selected from the data on how the transfer (or dis- 
tortion) function behaves in the B, region. In some cases these data are more readily derived than those for the 
frequency characteristics. 

2. The Necessary Conditions for Precise Reproduction 


1. Reproduction structure. We derive the reproduction structure for a;/(t) when the transform A;,(p) and 
K(p) are rational functions of p having no poles at p= 0. 


In this case the transform of ;_(t) is given by 


Aout(P) = 


* This interpretation assumes that one can effect a Laplace transform for any of the probable inputs in the 
class, and also that any such transform is a rational function of p, with all its poles lying ina finite re- 
gion of the p plane. (k) 

** In defining Q it is assumed that at all discontinuities in a, (t) and in its derivatives Qin (t) there are two 
equally likely values: 


a”) (th = tim oft (t + At) 


k 
(= _lim aff) (t — 


| 


and the output & ,,,,(t) itself by the Riemann-Mellin formula 


a+joo 
Gout (¢) = out (p) eP'dp. * (10) 


Since Aoyr(p) has no special points other than the poles lying in a finite area of the P plane, (10) can be 
represented as a sum of residues for the poles of Aju +(p)eP*: 


+2 lim [Ain (p) K (p) (p + + 


dp 
iy! i opt 
) [An (p) K (p) (p + 15)" 


Here a; is a pole of Ajn(p) of order r;, 6), is a pole of K(p) of order s x’ ¥j is @ pole common to A;_(p) and 
K(p) of order qj (for z values i and k; qj = i+ $k), Ny is the number of sels. of Ain(P), Nz those of K(p) and z 
the number of poles common to Ajp(p) and K(p). 


Comparing (11) with how precise servo and recording systems reproduce a; (t), we have 


inl 


K (0) r 
out (t) py il in (P) (Pp + %) te 


It is readily showf that the required correspondence between Oo,;(t) and &jp(t) can only be provided in a 
real system via the first term in(11). This term will henceforth be termed the forced (constrained) output com- 
ponent and denoted by @,,,,,(t)c. 


The second term in (11) describes the free motion of the system, and is hence denoted by Oj )_(t)- 


Finally the third term in (11) expresses the complex resonance state™ to which the system goes overon ap- 
plying oj, (t) to the input; it will be termed the resonance component, and denoted by &oyt(t) res. 


Hence 
our (t) = + “our + (res: (12) 


2. Reproduction error, The reproduction properties are evaluated in terms of the deviation between the re- 
productions of a;,,(t) by the real system and by a precise reproducing system. 


For a servo the deviation is 


Ao (t) = %our (t) — % out (4) 
which we term the instantaneous followup error, while for a recording system 
Ax (t)= (2) 


* In this case o must be greater than the real part of any of the poles of Aoyt(P). 
** Strelkov (2, 3] first used the term ‘complex resonance’. The ordinary frequency resonance is a particular case. 


N, 
K (0) 
2 dp 
| 


is the instantaneous recording error. 


Hence from (12) 


Ao (t) = (4) — %our + + out (Ares, (13a) 
Ay (t) = (Qe out our (Of + Mout (Ores: (13b) 


The first term in (13a) is the forced component of the followup error 4p (t),, the first in (13b) the forced 
component of the recording error A, (t),,. 


(13a) and (13b) imply that the free and resonance components of the followup and recording errors are identi- 
cal, i.e. 


A(t)¢ = Ao(t)p = Ar(t) 
A(t) reg = Ao (res = Ae (tres = out (res 


It is readily shown from the structure of 49(t) and that the components Au and 
Ayr)( res differ markedly from one another in orgin and in time variation. ; 
It is thus desirable to examine each error component (followup or recording) separately. 


As metrological criteria for evaluating the reproduction accuracy for a given class of inputs we can use the 
maximum values of the Tchebycheff criteria in the class: 


Kory = sup |Age)(t), = max (14a) 

Ar |, i. e. Ag =maxAr, (14b) 

Ares= sup |Agey(t)e |, max Ayes, 


and also the minimum damping rate = 


3. Subdivision of the conditions for precise reproduction. The nature of the correspondence between the 
class of inputs and the system suitable for reproducing them determine the conditions for precise reproduction. 


We must here differentiate the necessary from the sufficient conditions for precise reproduction. 
The necessary establish that the system is in principle suitable for the given class of inputs. 


The sufficient establish that the system is suitable for reproducing the given class of inputs with a set ac~- 
curacy. 


Two types of condition (precise followup and precise recording) can be differentiated in terms of the type 
of correspondence required between input and output. 


* From the practical viewpoint two factors can with advantage be used to specify the free motion: 
1, The Tchebycheff criterion eo | Gout (t)¢|, being the maximum value of the modulus to the 
oo 
input O&jp(t)= M1(t) in the interval 0 =t =o. 


2, The response time of the free motion, Tf, i.e. the time during which | Gour(t)f| does not exceed some 
preset value K(0)Meés 


KO) Mig = (Te = SUP (Op: 


4. The first of the necessary conditions for precise reproduction. (11)can be used to show that simple in- 
put reproduction is possible only in a stable reproducing system.* We can only speak of more or less precise re- 


production for a given class of inputs when the rate at which the free motion dies away substantially exceeds the 
rate at which any of the inputs does. This is the first of the necessary conditions for precise re production (follow - 
up or recording) for a given class of inputs. If the line Re(p) = —a— is the right boundary of the region of special 
points of K (p) (region P, in the figure) and Re(p) = — b — is the left boundary to the region Pq,the first of the ne- 
cessary conditions for gocies reproduction can be put as 


—a<—b. (15) 


If this is fulfilled in region Py there are no poles in the transfer system (z = 0), so no complex resonance 
occurs with any of the probable inputs in the class. Thus when the first of the necessary conditions for precise re- 
production is fulfilled 


Ares = max | A (tres | = 0, 
Moe) (4) = Airy 


5. The second necessary condition. Mere fulfillment of the first necessary condition does not ensure that 
the system is in principle entirely suitable for reproducing the given class of inputs. This only occurs if, together 
with (15) being fulfilled, the Tchebycheff criteria 


= sup |, Ar = sup {A;(t), 
are suffiently small for the given class of inputs. 
Since when z = 0 


4 
= i pt (a) 
Ay (t) c. Py dp’ (p) Ain (p) (p + aj) 
_KO rep! 
Ar (te = (p) (P) (Pp + (b) 


it is readily seen that A(t), and A(t), will be sufficiently small only if the moduli of « 9(p) and x, (p) 
and of their first (1; = 1) derivatives in region B,, are sufficiently small. 


This is the second necessary condition for precise reproduction, 

Let us delineate region Py, in some regions Hi, and A in the distortion planes x= o(p) andx, =K,(P), 
respectively (see Figure). 

Formulas (a) and (b) above show that the minimum 4), and A, correspond to x 9(p) and x , (p) be- 


ing uniformly continuous in some closed region containing region P,,, and t to there being no poles in the transfer 
function in the region P,- 


This in turn means that a necessary condition is that the transform of the | * region onto the appropriate 
distortion planes should be compact. The less the radius Ho of the circle encouiig H, the less Ay, for the 
given class of inputs; the less the radius y, of the circle enclsoing Hy the less A , .. 


* Korn first discovered this as a required element in a problem of his on integral synchronization, i.e., a prob- 
lem,of precise reproduction. Kharkevich [5] has given a fairly comprehensive review of the papers on distortion 
in reproduction. 


| 
« 


3. The First Sufficient Condition for Precise Reproduction - Precise Followup 


1. Structure of the forced component in the output. We assume the class of inputs and the reproducing 
system to comply with the necessary conditions. Then z = 0; hence 


We also suppose that the B,, region lies in the uniform convergence zone of the series 


Or, an 


k=0 


which represents che transfer function for the system within the circle 


4 
lpPI<R= 
Tim 


A uniformlyconverging series can be differentiated an arbitrary number of times. Hence we can substitute 
(17) into (16): 


p+ dp 


Supposing it permissible to change the order of summation, we have 


out (¢) kl K (0) x (18) 
ri op! 
x lim — [P*4in(P) (Pp + %) te”). 
p-—ai dp 


| J 
WD 
| 


The meaning of 


— apt [p*A,.,(p) (p + 


N, 


[Ain (P) (Pp + Se? = (19) 


and that Ajn(p) (p+ &j) i opt is continuous for the ares variables P and t in the region of some point aj 
for the interval 0 =t < oo. (19) implies that 


(t) = lim [A in (Pp) (Pp + 


or, in accordance with Schwarz 's theorem [10],on shifting the differential operator 


N, 


im1 


The resulting expression enables us to put (18) into the form used in servo systems theory [6, 7]: 
k=0 
| We shall prove that this formula is correct for a class of inputs more extensive than that used in its deriva- 
tion. 
Suppose that we apply to the input of a system with the transfer function of (3) an arbitrary bounded input 
Bi, (t) which has all its derivatives finite for 0 =t =o. If the initial conditions are zero the forced component 
in the output 6 ,,(t), is a particular solution to 
Bout (He (te +++ (te = 
= K (0) B + K (0) ampi™ (t). 


This solution is sought in the form of a series 


Pour (t) = 2) Cx8™ 


which can be differentiated n times. 
Substituting this series into the differential equation shows that when 


Cr = K™ (0) 


it becomes an identity.* 


* To prove this it is best to use (24a) and (24b) below. 


is revealed if we recollect that 
| 
28 


This proves that (20) is correct for classes of inputs of arbitrary configuration in the B, region if the neces- 
sary conditions for precise reproduction are fulfilled.* 


2. The first sufficient condition for precise recording. We stated above that in evaluating the forced error 


component for a given class of inputs we should use the parameter Ay, = Max sup | Ao(tie | , the metrolo- 


gical analog of which is the absolute error of measurement. It is also appropriate to use parameter égo = 
the metrological analog of which is a relative error of measurement. 


We now evaluate 4p, and € 9, for the case where the class of inputs is defined by Q and M. 
By definition Ao (t), =&oyr(t)e —%o our (¢); hence, from (20) and (1), we can put 


Ben, 
K(0)M’" 


bo (de = 3) (0) 


k=1 


Substituting the resulting expression into (14a) we have 
k) 
4 ain () 
=k (0) M max sup |. 0) 


Evaluating we get 


(k) | 
Acc <K OMS | (0) |max sup | 


(21) 
(k) 
Comparing max sup | ae | with (9) we see that 
(k) 
max sup < at. 
This enables us to put (21) as 
be <K(0)M > 0) 10%, 
koi 
from which** 
<>) (0) | (22) 


k=l 


* It can also be proved that (20) is correct for linear systems with distributed constants. 
**It can be shown that, when not less than the first m derivatives in the given class of inputs are defined, we can 
use to estimate ¢€,,. an inequality 


i 


—jeo 
Equation (22) is clearly a special case of (22a), since m = oo and jim Lom = = lim, | "| 


Certain difficulties are encountered inusing (22); we have to compute the series 


(0) 1 (23) 


Hence, practically speaking, in determining €,, it is best to use inequalities such as €9, = nj, where nj; is one 
of the calculable majorants of (23). 


The majorants {n} enable one to express the first sufficient condition for precise reproduction in the form 
Ni = €T,, where €fc is some preset value. . 


Compliance with this inequality is sufficient to ensure that the relative deviation of & oy (t), from &oy_(t) 
for any of the inputs in the given class will not during time 0 <t = w exceed €j¢. 


Strictly speaking, we can construct an infinite number of calculable majorants to (23), and hence an infinite 
number of forms for the first sufficient condition. But in practice only those majorants are of interest which have 
sufficiently simple formulas when (23) is evaluated sufficiently accurately. 


3. The majorants of > 4 | x(*) (0)| Q* . We will derive the majorants for certain cases of determining 
xp) and K(p). 
Let the distortion function be 


_ —by) p+ (a) — by) + ... 


defined by the values 


It can be shown [12] that the coefficients in the Maclaurin series of the distortion function B, = een (0) 
are related by recurrence formulas 


B, = AiAy_, 4+ AeA, (24a) 


where 


for i<0, 


a,=0 for i> m, b,=0 for 


(24a) enables us to represent (23) as an infinite sum ) | B, | Q* , each term in which can be evaluated 
from 
| Q= | | Q, 
| By | | A; || Ar | +} | 0°, 


We hk 
k=1 i—1 k=1 


a,=0 for idm, b,=0 for i>n). (25) 


| 
| 
| 
| 
pBy | Q* Ay] | Ay_, | 1 As | | OX +... A, 0". 
| 
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ted 


Remembering that when ay= a,=...= a), = 0,B, = Ay we can, by using (24b), put 


k=1 k=1 
from which 


1410" < 


i=—1 


Substituting the estimate of | 4, | Q* into (25) we get 
k=1 


| 


(26) 
|b, al 


Then the €9, is given the best for the class of inputs defined by Q and M, and for the set of reproducing 
systems defined by the values | by], ...,| —byl. | ag —agl 


Similarly, it is easily shown that 


06 
|b, | at 
(—1)! | 4,1 
% (p) = 


Analysis of 1) shows that it gives (as well as €9¢ ) an estimate of the maximum values of | « 9(p)| 
in the circle | p| = Q: for the given class of inputs: 


= 
i— 
This implies that 
€oc (27) 


* There is every reason to suppose that €9¢ © Hg retains its validity when x ¢p) belongs to a class of funct/ ons 
more extensive thanthe rational ones used here. 


In particular, we can prove it correct when x 9(p) is a holomorphic function of a single sheei within the 
circle | which contains the By region. 
lim | 
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= 


Estimates of (23) of forms other than that of (26) can be made using combined parameters for the re produc- 
ing system. 


Such parameters are® 


Ry = 


1 
k) 


(k) 
D, -~ Xu (0) |, 
kt 0) | t) 


1 
We now derive the majorants of interest to us by using the grouped parameters enumerated above. 
a) Let the distortion function be defined by Rj, and Ra, and by the numbers m and n. 


By definition, Rj, is the distance between p= 0 and the nearest zero of b(p), while R, is the distance between 
p= 0 and the nearest zero of [A(p)+ 1] = &(p). 


It is readily shown that the coefficients in b(p) and 6 (p) are relatedto Rp and Ra by inequalities of the type 


1 i 
b 


Consequently, 
® ai n 
Q n 
314,10! <3 -1. 
int imi 
Substituting these in (26) we get 
“oc (28) 
* The parameters Rp and Ra, to determine which Cauchy-Adamar fomula is used [9], are the convergence radii 
A(p) b (p) b (p) if the 
of the Maclaurin series for x 9(p) b(P) and A(p)+ 1 8 (p) » Tespectively. It can be shown that 


coefficients inthe polynomials b (p) and & (p) are everywhere positive, it is permissible to replace Rp and Ra in 
all the estimates of (23) given below by the lower bounds to the zeros of b(p) and §(p),respectively. From a 
certain theorem [11] the bounds are 


| 
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b) Assume a distortion function defined by D,(a parameter). 
By definition, | B, | =< Des hence 


kewl 


If the distortion function is defined by the coefficients in a(p) and b(p) the By needed to determine D, can 
be found from 


4; A, As A, 
0 
= (—1)'—"Dy, (iS k #0), 
by be 
0 01 2B 


where 4; = 0 when i> I, bj = 0 wheni> n, aj = 0 wheni> m. 
c) Letting the distortion function be defined by Q, we can write 


| (0) |< 


k=1 


(30) 


Q» may also be found from 


d) Let the transfer function be given via the parameters L and Rp, and also by the numbers of poles by orders. 
By definition 
eee eee 
K(p) = 1K 0)| + + 
where sup Ix |; = 1, N is the total number of polesof K(p) and 8; is the order for pole 8; . 
Sy 


Hence, here, for k = 1, we haves 


4 


| 
+ Chas (+--+ )+ 


1 N 


| 
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i i 


44 N 
% N-1 i 
(fiat 


| (where Nj is the number of poles of order not less than i) we can put 
N N N, 
Substituting this estimate of | B, | into (23) we get 


This estimate of can be somewhat simplified. Nj/ Ny Ishence 
siss 


| LN (R,—Q)’—1f (R,—O)* (R,—Q)—1 
| 


In the special cases = 1 we have 
Q LN 
“ce 


Analysis of the estimates of €9_ shows that they all have certain common properties. 
1. All are absolutely precise for the input» 


a, (t)= Mi(t) (Q=0). 
2. All are absolutely precise for a precise servo.* 


3. Like (26), (28)-(31) are the best for a given class of inputs and for the set of systems defined by the grouped 
parameters used in the estimates.** 


These properties indicate that our estimates of €9, above are compratively of high accuracy. The accuracy 


is the greater the smaller Q, the smaller the size of the P,, region and the more monotonic the distortion function 
in this region. 


| A further property of the estimates of €9¢ is to be noted. 
Considering the majorants of (23) for values of Q which ensure et oe is complied with, and re- 
membering that in most practical instances 0 Land | (1 <k =<), it is readily seen 


that for any of our marjorants 9; (9) lim [nj (Q)). 


* For a precise servo ay = a,=.. . = am, = 0, =0 or aj= bj (i= 1, 2,...,m), Rh = RA= 
= and Dp = Q,=L* 0. 


** Thus when reproducing ojn(t) = Me~* (a class of inputs specified by M and ) with a system of distortion 


a function x ¢p) ——— (a set of reproducing systems specified by Dy) we have sup | (t) a = 
= - Ost=@ 


Dy = 
1 - QD, ns 


| Using the fact that 
| 
| 
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Thus to perform a preliminary check on the correspondence between a given class of inputs and the reproduc- 
ing system we may use the relationships 


toc 
Of these, one of special interest is * 
from which €9, can be estimated from the experimental phase-frequency response. 


4. The First Sufficient Condition for Precise Reproduction~-Precise Recording. 


1. Structure of A. (t)c We have already shown that the forced components of the instantaneous record- 
ing error can be calculated from 


A, (the = four (t)c — (t). (33) 


Provided the necessary conditions for precise reproduction are complied with, the first term @,,,,(t), may 
be shown as the functional series 


x” (0) (t) our(te (34) 


k=0 


while the second, Op oy_(t) (0) (.. 
0 
K(0) (t) = ult) 
Substituting (34) and (35) into (33) we get 


= > (0) = dere (P(t), 
k=0 


for k<i 
= |e K0) -| 


for k>1. 


“K@) K(0) 


* There is every reason to assume that * 9(p) retains its validity when k belongs to a class of functions more ex- 
tensive than the rational ones used here. 


In particular, we can prove it correct when « 9(p) is a holomorphic function of a single sheet within the circle 
which contains the P,, regions 


lim 


1<k- @ 


R 


The proof of this involves using the properties of single-sheet functions which are holomorphic and normalized 
within the circle | z 


2. The first sufficient conditon for precise recording. By using a proof similar to that given for (22) we can 
deduce the following estimate for the limiting value of 4, ¢ for the given class of inputs: 


k=2 


A 
< Di Cea". (36) 


The practical use of (36), just as with (22), involves certain difficulties in calculating the series 


| Q*. (37) 


These difficulties are removed by replacing series (37) appearing in (36) by some one of its calculable ma- 
jorants. Let €; be one of these. Then the first sufficient condition for precise recording can be put as €9¢ = 
c, where €, is some preset number. 


Compliance with this is sufficient to ensure that the system will reproduce any of the probable inputs defined 


by M and Q which has values of Bren sup | Ar(t)c| , which do not exceed K (0) Mé» c. 
stso 


3. The majorants of >) |C, |". By definition, the distortion function 
k= 


is holomorphic within the circle | p| < R. Hence within this circlex,(p) can be represented as a Maclaurin series: 


k=0 


k=2 


The coefficients in this are functionally related to those in (23). In fact, 


Using this, together with (24a) and (24b) we can derive recurrence relations for the set of coefficients Vt: 


Ay = — — — (39) 
A, => a; am bi, 


where 
0 for i>0, a,=0 for i>m, 
for i=0, b=0 for i>dn. 


(38) enables us to represent (37) as an infinite sum with terms defined by inequalities of the type 


k=2 


+ (1+ 314,00"). 


kel 
Herel , as before, is the larger of m and n. 
Remembering that 
n 
> a8 
>) 
we can put €_ ¢ in the form 
2° 
i— 
imi 


We can use other majorants than the , to estimate €, ¢, i.e. the ones constructed by using the grouped 
parameters D,, Q,, etc., of Section 3, 


The majorants, in particular, are 
& =n, 
“TOD, 


& 


They are constructed on the same principles as the 9;. 
Certain typical features of the estimates for¢, are worthy of note. 


1, For all majorants lim = 0, while lim my 0, Hence for @ sufficlently small 
ti 


2. For all majorants }g ¢ lim §; = ;. Hence for small values of 


» 
| = =| 09) = |a, 
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This then gives | 


the requirements as to precise followup and recording practically coincide. 


These properties of the estimates for the ¢, , enable us to set practial bounds to the utility of this parameter 
for evaluating the reproduction propercies of any particular system. 


5. The Second Sufficient Condition for Precise Reproduction 


1, Structure of our (t)f. We will derive the structure of the free component for the case where the transfer 
function and the class of inputs comply with the necessary conditions for precise reproduction. As before, we sup- 
pose that the transfer function and the transforms of all probable inputs are rational functions of p, all the poles of 
which lie in a finite region of the p plane. 


Here z = 0, so from (11) we have 


N 
= 2 lim = [Au (p) K (p)(p+ By)**er*]. (41) 


p—>— 


Apply Leibnitz's formula to the components in the RHS in (41) and, passing to the limit, we can put 


= De™ Dart! Fr * 
k=1 j-0 


a 


x lim K Pay 
Since, in accordance with the notation previously adopted, 

1 


= K (0) Li, 


8-1 


N 
4 
Sour (t) ¢ = K(0)L de Bx X 
kel 


When is reproduced the component o,,,,,(t), occurs whenever (t) or it derivatives aK) (4) show 
discontinuities. 


In the general case an input with discontinuities at ty, tz, . . ., ty cam be represented as a sum of inputs: 


in(¢) = Gin (t tx), 


where each component, plus: its derivatives, is continuous for t, = t = o and precisely zero for t < ty. 


In this case the free motion will occur at t = t,, so from (42) we have in each case? 
when t < ty 


As = ous =0: 


a 
wi 
He 
(42) 
it. 
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when t 2 ty 


N 


k= 

t,) > Pi. Ain (— Bx). 
i-0 


2. The second sufficient condition for precise reproduction. At the instants when | our (t) ¢ is compar- 
able with |@our(t)| and | aout (t)| the system does not, strictly speaking, fulfill its function. 


This is because we have required that the free motion be rapidly damped, which is imposed because of the 
necessary conditions for precise reproduction. 


It is best to use Tf to evaluate the damping rate; it is given by 
K(0)Mt; = su t 
If Tf = Tg and the first sufficient condition are complied with the system reproduces the set of inputs very 


well. In this case the metrological parameters €9¢ ,€rc¢ » 4nd Tr will not exceed some preset values €oc 
€rc + Tftespectively. 


The second sufficient condition is thus that the inequality Tr = T, be complied with. 


3. The majorants of &out(t)¢. Using (42) with Ajn(p) = Ain(p) = 3 we get 


-1 
N 


-1 i=0 
Hence, 
<t<e 2 Pe 2 2 


As this expression is complicated it is very difficult to determine Tg. Let us consider ways of simplifying 
it. 


Remembering that 
| Be | he 
from (42), write 
RE —1 1 RR 
Re—1 RS = F(t). 


The majorant F,(t) is correct for any arrangement of the N zeros of b(p) of degree n = Nys in the semiplane 
Re(p) <—a (excluding the area within the circle | p| = Rp). 


_ 


This favorably differentiates F(t) from F(t) ; 


proposed by Fel'dbaum for certain classes of distribution for the zeros of b(p) of degree n in the semiplane Re(p)=< 
ss ~a (5). 
However, the majorants of greatest practical interest are 


1 


= = —(a—1)t 
F,(t) K (0) MF ze for a>i 
and 
F, = K LMN 4 — 


= K(0)MF,t-' for a>0, 
which are suitable for finding Ts rapidly. 
The first is constructed using 


4 


which is correct for Rp > 1 andt> 0; the second by using 


i+1 i+1 


thee [=] ‘for >0 


into F,(t) we also get an estimate Mg = our 
— 1) + Ri — 1 (s—i)! (2) 


Hence F,(t), F(t) and F(t) enable us to show the second sufficient condition either as 


1 PF; 


We finally observe that F,(t), F,(t), and F3(t) remain the majorants | a5,,(t)¢| when Rp is replaced by a. 
If the coefficients in b(p) are always positive we can also replace Rj, by the lower bounds to the zeros of this poly - 
nomial ¢ 


1 
I 
1 
t 
| 
vi 
: 
for t>0. 
Substituting 
or as 
Ps 
a>0. 
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SUMMARY 


1. The methods given for evaluating the reproducing properties of linear systems enable one to decide on 
the correspondence between the system and the class of inputs which will ensure a given precision (in followup 
or recording). ‘ 


2, These methods may be of value in solving problems of analysis, synthesis and correction for measure- 
ment, servo and control systems, ~ 
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THE PROBLEM OF SYNTHESIZING LINEAR DYNAMIC SYSTEMS WITH 
VARIABLE PARAMETERS 


A. M. Batkov 
( Moscow) 


A method is presented for determining the differential equation of a linear 
dynamic system with variable parameters, when the impulse transient response is 
given. 


One of the important problems in the synthesis of automatic control systems is the determination of a de- 


sign to implement a known impulse transient response. 


In the case of systems whose processes are described by linear differential equations with constant coeffici- 
ents, this problem is solved by determining the transfer function corresponding to the known impulse transient re - 
sponse of the system, and then representing this transfer function in the form of a product of elementary factors. 


If the system is described by a linear differential equation with variable coefficients, then its transfer func - 


tion depends not only on frequency, but also on time [1], and this, in the majority of cases, renders it difficult 


to determine the design for implementing the system. Because of this, one usually proposes to determine the de - 
sign of systems with variable parameters by starting out from the differential equation characterizing the system. 
Below, we consider the problem of determining this differential equation in the cases where the impulse transient 
response is given.* The questions involved in determining the design corresponding to the equation are considered 


in the examples. 


Let there be given an impulse transient response W(t, 1) of a linear dynamic system with variable para- 


meters. We assume that the differential equation for this system has the form 


L(p, t) u(t) = M (p, t)o(t), 
where the operators L(p, t) and M(p, t) are defined by the expressions 


t=0 


and p= d/ dt, u(t) is the output and v(t) the arbitrary input variables of the system, the variable coefficients, aj (t) 


and bj (t), are continuous functions of time, and n and m are unknown constant integers. 
One of the possible methods of determining these quantities is expounded below. 
It is known [2] that the output variable of System (1) may be written in the form 


u(t) = (W(t, dt, 
0 


* This happens, for example, in problems of synthesizing optimum filters. 
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(1) 


(2) 


(3) 


i 

= 


where W(t, r) is the impulse transient response of the system, and is assumed to be known. 
From Expression (3) it is easy to determine the initial value of u(t) and its derivatives; 
u(t) l, ~= 0, 
u’ (t)|,_, = w(0) (0), 
W(t) = 2w’ (0) (0) + w(0)v’ (0), 
(t) = kw(t— (0) v (0) + wht—2) (0) v’ (0) +--+ 
+ kw’ (0) (0) + w (0) (0). 
The problem consists of determining, in the form of (3), the solution of differential Equation (1), correspond- 
ing to the initial Conditions (4). 
Writing the general solution of the homogeneous equation 


L(p, t)u(t) =0 (5) 
in the form 
w (t) = cy (t) + (t) (t) +++ ++ (t) (t) (6) 
and applying the method of varying the arbitrary constants [3], we obtain a non-homogencous system of n equations 
in the derivatives, cj(t),..., cy 
+ tat. + =0, 
+ +... + = 0, 


+ +... Mey M(p, t)v (t). 


Since the determinant A(t) of this system is the Wronskian of Equation (5), it is different from zero and, 
consequently, System (7) has a unique solution, which may be written in the form 


(7) 


A 
=) t)v(t)dt+%. (8) 
(t) 
The constants y; are determined from the initial Conditions (4). In order to so determine them, we sub- 
stitute u(t) and its first n-1 derivatives from (6) into (4). Since the non-homogeneous system of equations thus ob- 
tained has its determinant equal to A(0), it, as with (7), has a unique solution 


u, (0).. Ow uy, (0) 


Since, from (8), 


* 


— 


then, comparing (8) and (9), we obtain finally, 


t | 


Consequently, a particular solution of the linear, non-homogeneous differential Equation (2) with variable 
coefficients and initial Conditions (4) may be written in the form 


t 
ain A; (0) 


where (t), + » Up (t) are particular solutions of the homogeneous Eqautions (5). 


Equating u(t) from (3) and (11), we obtain relationships for determining the particular solutions, uy, ug, . . . | 
. +» Up, and the right member, M(p, r), of the desired differential equation in the form ‘ 


t 
ain A( (0) 


Since the systems considered in this article are described by linear differential equations with variable co- 
efficients then, as a consequence of the theorem as to the existence of solutions of these equations [3], their im- 
pulse transient responses may always be given in the form * 


W (t, t) = D) g(t) u(t), (13) 


i=1 


where the g(t) are linearly independent functions, constituting a basis for the homogeneous differential Equations 
(5), corresponding to (13), and n is the order of the equation. 


If, in (12), we make the substitution, g(t) = u(t) (i= 1,2,...,M), we obtain the basis for the homogen- 
eous differential equation, which coincides,to within a factor independent of the derivatives,with the left member 
of the desired equation. 


Employing well-known methods [3], we determine the left member of the equation by expanding the (n+ 1)- 
-order determinant: 


Dit)e 


wim 


Oo 


With Condition (13), and as a consequence of the linear independence of the functions uj(t) (i= 1,2,....9), 
the relationships (12) can be written as a system of n equations of the form 


t A* (0) 


from which the form of the left member of the equation, i.e., M(p, t), is determined. 


* We do not consider here the questions regarding transforming the given function W(t, r) to the form (13), and 
the corresponding identity of the impulse transient response of the system and its differential equation. 
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The computations are significantly simplified if it is necessary to determine a system of n linear differential 
equations of the first order, with variable coefficients, corresponding to the given impulse transient response of 
Form (13). 


In fact, taking W(t, r)as the sum of n functions, Wj(t, r) = ;(t)p4(r),in accordance with the reasoning 
above we obtain a system of n first-order differential equations, the operators of which are defined by expressions 


d 


= (i = 1, 


where M,* (p, t) is the operator conjugate to M, (p, t). 


It is assumed that these equations can be physically realized, i.e., that the system with the impulse transient 
response (13) can be implemented by n linear, first-order dynamic systems, connected in parallel, This is a very 
convenient approach from the point of view of applying linear analogue devices, 


The proposed method is completely applicable to systems with constant parameters. The method perm‘ts 
one to determine the form of the differential equation and, consequently, the design of the system in the casrs 
when the functions entering into W(t, r) are given in general form. It is clear that in such cases, even when the 
system has constant parameters, the method of determining the structure and parameters of the system from the 
transfer function is inapplicable. An example corresponding to this case is considered below. 


Simple examples to illustrate the application of the method. 
Example 1. In [1] there is considered a system, the impulse transient response of which is equal to 


W (t, 2) = (1 + exp {—[¢— 2) + (sin —sino)]} 


In this case, n= 1, and g(t) = exp}- (t+ sinwot ) t= A(t). 
The determinant corresponding to (14) will, in this case, have the form: 
Sane} 


~exp{—(#+ sin agt)} (1 + 9.008 


Expanding the determinant, and discarding the multiplier exp \- (so Le sin wet)f, we obtain the left 
member of the desired differential equation: 


L(p, t) u(t) =u’ + (1 + p cos wot) u. 


For the determination of the right member, we find 4(0) = 1, Af (0) = 0. Equation (15) then takes the form 
t 
cxp £ sin age)} (p, 1) + age) = 0. 


M (p, t) = 1 + pcos aot. 
Thus, we obtain the equation for the system: 


+ (1 + 008 u (1 + 9.608 v(t). 


From this it is clear that we have a first-order system which can be realized by an RC filter, with variable 
resistance 


i 


As determined in [1}, the approximate transfer function for p « 1 will be 


jo 


je 6 


Whence the imposibility of drawing any conclusion as to the design of the system. 


Example 2. ‘It was recently shown [4] that for optimum filtering of a random signal in the presence of 
white noise, the linear system with variable parameters would have an impulse transient response of the form 


t 


W(t, 
N* + do 
0 


where g(t) is a known function . 


We shall determine, by the method expounded in this article, system design and parameters corresponding 
to W(t, r). 


Obviously, 


n= fandp, = 
N? + do 
0 


t 
+ 20] — 


This gives us the left membet of the equation 


t 
L(p, t)u(t) =u’ — : 


(m+ 


For determining M(p, t), we employ (15), obtaining the equation 


@*(0) do 


M (p, 


an 
| 
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M (p, = 
do 


Consequently, the system under consideration is described by the differential equation 


fe 


and may thus be realized, for example, by a first-order RC filter, in which R and C vary with time according to 


t 
N* + § g*(6) do 


=, 


Example 3. We consider now an application of our method to the determination of the differential equation 
for a system with constant parameters, 


Let 
W (t, t) = cos way (¢ — 1). (16) 


Transforming cos w»(t= by well-known trigonometric formulas, we obtain n= 2, g(t)= COS Wot, 
x(t) = e 8" sin w,t. Setting up the determinant D(t) according to Formula (14), we find that 


L(p, t) = (p+ a)? + 


Since A(t) = 


we get the following system of equations for determining the right member of the equation; 


£08 wytM (p, 2) = 0. 


Integrating by parts and adding, we find 


t 
gin (¢ — 1) [M (p, — 0’ —av(x)] dr 


Thus M(p, t)= p+ a and the desired differential equation has the form: 


It is easy to verify that the direct application of the Laplace transform to the impulse transient response, 
(16), leads to an analogous result. 
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TELEMETERING SYSTEMS WITH PULSE=-CODE MODULATION 


G. V. Burdenkoyv 
(Moscow) 


High-speed pulse -code telemetering devices are investigated. The possi- 
bility is demonstrated of constructing telemetering deivces with magnetic elements 
having rectangular hysteresis loops, together with crystal triodes and diodes, The 
basic parameters of such devices are cited, and an estimate is given of the accuracy 
of telemetering. 


Experiments in the use of telemetering devices in power systems show that the unstable functioning of these 
telemetering devices is to a large measure conditioned by the influence of interference in the channel, interfer- 
ence which is particularly noticeable when the channel is complex or extensive. Hence the significance of ques- 
tions of stability under interference of telemetering systems, and of questions of the independence of telemetered 
findings from variations in the parameters of the channel. 


As is well known, among the telemetering systems most stable under interference is the pulse-code system, 
the functioning of which is the most independent of the channel's state, since the discrete method of telemeter- 


ing does not entail an accurate reproduction of pulse form, but requires only a determination of the presence or 
absence of a pulse (see Fig. 1). 


Signal 


‘Time 
Fig. 1. Recovery of the code on the receiving end. 
I) Clipped below; II) Clipped above. 
Shadings mark the recovered pulses. 


The very simplest principle of code formation is that of numerical coding, i.e., the construction of a sequence 
of combinations using one of the standard number systems, e.g., decimal or binary. Binary coding is usually em- 
ployed in high-speed pulse-code sysvems, since the binary code, for the same accuracy in measurement, advan~- 
tageously distinguishes itself from the decimal code by significant decrease in the necessary number of pulses, thus 
permitting a desirable decrease in the magnitude of error in discreteness for sufficiently high speeds of function~- 
ing. The complexity of implementing the coding and decoding must be mentioned as a disadvantage of these 
methods. The most well-known forms of such coding are: 


a) Pulse-count coding; 

b) Inverted pulse -count coding; 

c) Coding by the principle of comparison with reading pulses; 

d) Coding by means of variations in the accumulator voltage; 

e) Coding using a comparison scheme with a multi-step comparison voltage [1 - 4]. 
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Transmission may thus be carried out either by qualitatively different pulses, or by the presence or absence 
of homogeneous pulses. 


Not all known forms of coding satisfy the requirements stemming from the conditions under which telemeter- 
ing must be employed in power systems (as to accuracy, stability, etc.). 


From this point of view, there is a good deal of interest in methods of coding using comparison schemes with 
multi-step comparison voltages (for brevity, these will be called “methods of multi-step equalization"), which 
were taken as the basis of development of pulse-code telemetering devices. 


Particular attention was paid to the choice of elements for the apparatus, including the use of magnetic ele- 
ments with rectangular hysteresis loops in conjunction with crystal triodes and diodes, with respect to the varia- 
tion of the element parameters which, in the given case, has little influence on the accuracy of the telemetering, 
since we are dealing here with bi-stable devices, analogous to relays. From this was envisioned the possiblity of 
modular devices, the number of modules, or cells, to depend on the required number of elements in the code, as 
determined from the inequality, 5p= 4(2(2" - 1))7, where §p is the discreteness error and n is the number of 
elements (pulses) in the code. 


The frequency F with which code pulses may succeed one another is a function of the channel band-width, 


Af. 


For our method of telemetering coding, as already mentioned, it suffices merely to distinguish the presence 
or absence of a pulse. Therefore, it is possible to limit the band-width, Af = 2F. 


Code transmission time (time for one cycle) is found from the formula 


ta" n(t, + + tg+ ti» 
where n is the number of elements in the code, t, is the length of one code pulse, tp is the length of the pause 
between code pulses, t, is the length of the pause between code groups, and t,,, is the length of the synchronizing 
pulse (the marker). 


In the device developed at the Central Laboratories for Scientific Research in Electrical Engineering (TsNIEL), 
the choice of parameters was; n= 6, F = 40 c/sec t, = th= tg tm = 3ty. The time for one cycle, T, = 0.2 sec, 


Fig. 2 is a block schematic of the transmitting device. The design has the following component parts: MV, 
a multivibrator, generating pulses for controlling code pulse distribution; EF1, EF2 and EF3, emitter followers for 
load-matching; D is atwo-way distributor of the pulses fed to the magnetic elements with rectangular hysteresis 
loops, I, I", etc., these pulses switching the coding cells in the proper sequence; T is a trigger assembly, for switch- 
ing the conductor assembly, comprised of crystal triodes; CA is the conductor assembly (in this circuit, electronic 
switching is performed by 6H2 tubes); ZD is the zeroing device, controlled by the sign of the difference between 
the measured (U,) and the equalizing (U,) voltages; CC is the conductor circuit; CB1 is the coincidence block, 
generating pulses for cutting off the corresponding trigger cell; CB2 is the block which isolates the output code 
pulses; SC1 and SC2 are shaping cascades for transforming the input pulses to rectangular pulses. 


Coding is accomplished in the following manner; 


As a preliminary step, the measured quantity, varying with time, is transformed to a rectified voltage, Uj. 
During the coding process, this voltage is compared with the equalizing voltage, U,. This latter voltage varies 
stepwise, the amplitude of the steps being proportional to the decreasing binary series, 2"~*, 28-2, 20-%,..., 2°, 
The zeroing device is controlled by the sign of the voltage difference, U, - U,. If this difference is positive, then 
the corresponding binary step is retained. For a negative difference, the corresponding binary step is not utilized 
in the formation of the equalizing voltage. In order to retain a positive difference while simultaneously switch- 
ing to the next binary step, a pulse is sent through the path of CB2, this pulse corresponding to the previous term 
of the code, 


Let us clarify the coding process by way of an example (Fig. 3). Assuming a six-valued code is being used, 
if the measured voltage U,= 45U, (Uy being taken as the relative unit), then the following steps must be included 
in the equalizing voltage: 2°Uy, 2°U, and 2°Uy. This will give the equalizing voltage U,= (32+ 8+ 4+ 1)U)* 
= 45U, = U,. Thus, code pulses 2°, 2 2", and 2° will be transmitted from time 2 through time 13 in the channel. 
Each cycle is preceded by a synchronizing pulse (marker) whose length is significantly different from that of the 


code pulses, The marker is produced by element M of the two-way distributor and by shaper SC2. Distributor 
elements CT and BP prepare the circuit for the following code cycle. One of the basic components of the trans- 
mitting device is the two-way magnetic distributor, working with elements having rectangular hysteresis loops [5]. 
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Fig. 3. Graph of the coding process. 


The comparison element in the transmitter design 
i Z is the zeroing device (Fig, 4). Its design is comprised 
of two magnetic elements, I and Il, with rectangular hys- 
teresis loops, and a trigger T with one stable state. 


4 +™% T to A current proportional to voltage U, flows through 
CB1 winding Wy and a current, from the conductor circuit, 

i {Wi proportional to compensating voltage U,, flows through 
winding Wy. Depending on the sign of the voltage dif - 
We fire ference Uy - Uge core I is or is not primed to be switched. 


If the compensating voltage is the larger of the two, 
then core I is primed and, when sensed via control wind- 
Fig. 4. Circui roing devi 
8 laos s ing Wc, Causes a pulse to be generated at winding W,, 
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which in turn transfers trigger T. This then provides a gating voltage to coincidence block CB1 and yields a pulse 
to switch off the memory trigger corresponding to the proper step of compensating (equalizing) voltage. 


Fig. 5. Circuit of the stepwise equalization voltage assembly. 


Core II plays a secondary part, serving to dissipate “parasitic” counter pulses between load winding Wh and 
anti-parasitic winding W,p. The core itself is prepared from type 65 NP material. 


It should be mentioned that ii 's preferable to use ferromagnetic material with low coercivity for the zero- 
ing device, thus increasing the sensitivity. 


The shaping cascades are made up of monostable triggers with type P1B crystal triodes. 


The assembly of stepwise equalization voltages is implemented by a circuit with electronic switches, the 
switching of which is accomplished in a definite order by the triggers corresponding to the code cells (Fig. 5). The 
electronic switch is normally cut off by the voltage taken off the shunting trigger triodes. 


_ The basic data for computing the magnitude of the resistance r_ are the currents I) and I, and the resist- 
ance Ry. Calculation of the common shunt is made by means of the formula, r= k(Rx + R,,,), where ly is the 
current flowing through the tubes, I, is the magnitude of the current through winding Wy of the zeroing device, 
Rx is the resistance of winding Wy, and 


keg = 2° + 224... 422, 


ry = ky (Ru t+ Rin), Te = ke (Ryt+ Rio) 
Tn = ken (Ryt Rio) — + 


Here, ko, . . . , ky are coefficients. 


The coincidence circuit (Fig. 6) of the transmitter is based on the principle of shunting diode switches. With 
such a circuit, one of the compared voltages feeds the input, and the other cuts off rectifier D,, shunting the out- 
put, The circuit is controlled by the coincidence of pulses taken from distributor element I and gating voltages 
from trigger T. With this voltage present across the voltage dividers R, R, (the rightmost triode of the trigger 
being cut off), diode D, does not pass current, and, with the load winding, no pulse voltage drop occurs across re~ 
sistor Rp ; the pulse is transmitted without loss to the pulse shaper via capacitor C. When the trigger reverses, the 
gating voltage is absent at diode D,, and a large portion of the voltage is dropped across resistor Rp (Rp >> Ry, )- 


The interaction of the sub-assemblies of the transmitter is shown in Fig. 7. 
Fig. 8 isa block schematic of the receiver. The design consists of the following blocks: A is an amplifier; 
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AC is an amplifier-clipper, used for two-sided clipping 
of the channel pulses; MV is a multivibrator which 
generates pulses to control the magnetic distributor; 
TS is the time selector, used to discriminate the longer 
synchronizing pulses; EF1, EF2 and EF3 are emitter 
followers; D is a one-way magnetic pulse distributor; 
CB is a coincidence block; T is the memory trigger 
block; CA is the conductor assembly; CC is the con- 
ductor circuit; OI is the output indicating unit. 


The one-way magnetic distributor employed in 
the receiver is composed of the same elements as the 

Fig. 6. Coincidence circuit for pulses from 

element I and trigger voltages. two-way distributor of the transmitter. The parameters 
for the magnetic elements are chosen to be the same 

as in the two-way distributor. 
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Fig. 7. Cycle diagram for transmitter functioning. 

A) Pulses controlling the distributor; B) Switching diagram of the 
distributor elements; C) Switching of the zeroing device; D) Switch- 
in diagram for the triggers; E) Pulses entering the shaping cascades; 
F) Pulses in the channel. 1) element is primed; 2) moment of 
switching; 0) element “void.” 
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Fig. 10. Interaction of memory elements K with the 
distributor elements. 


The pulse time selector (Fig. 9) consists of an intermediate amplifier IA, a trigger with one stable state, 
T, and emitter follower EF. The selection principle is based on the selection of a definite time constant at the 
input of the intermediate amplifier. Such a selector is frequently called an integrator. If a short pulse of length 
Ts and a long pulse of length ry appear alternately at the input to the selector, then at the output of the ampli- 
fier we obtain the ratio of the peak voltages, which is approximately equal to the ratio of the lengths. 


The excitation threshold of the trigger is so chosen that the trigger will fire only for pulses of length r; . 


The remaining assemblies of the receiver are duplicates of units in the transmitter and, hence, will not be 
described again. 


For decoding, the distributor D is controlled by a pulse generator which functions synchronously with the 
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generator on the transmitting end. Synchronization is implemented by the synchronizing pulses (markers) and by 
subsidiary pulses entering into the code combinations. 


Moreover, start-stop distributor-synchronizing pulses, obtained by use of the "void" distributor elements 


I-VI and the preparatory first element B, which discriminates the synchronizing pulse from the code pulses, are 
provided on the receiving end. 


Code and synchronizing pulses entering the channel amplified and clipped, impinge on one of the windings 
of magnetic memory element K and prime it every time. Sensing of element K is executed under control of 
the multivibrator simultaneously with the switching of the magnetic distributor 180° out of phase with the code 
pulses. Thus, after a synchronizing pulse impinges on the load winding of element K, a pulse is formed, simult- 
aneously with the formation of a pulse at element B. If the first code pulse to arrive were 2°, then in the next half 
period, pulses would be formed both at the load winding of K and at the winding of element 1, etc. The inter- 


action of these elements is shown in Fig. 10, Winding W; 9 of element K and all distributor elements I through VI 
would be switched in the reverse direction. 
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Fig. 11. Cycle diagram for receiver functioning. 
A) Code pulses; B) Control Pulses; C) Selector pulses; D) Switching 
of the memory elements; £) Rectifier switching; F) Trigger switching. 


Winding W, 4 is an element of coincidence block, CB. If a memory element were primed by a code pulse, 
then sensing would cause a pulse to be formed only at distributor element winding W, 1, but would suppress pulse 
formation at winding W, 9. Thus, a pulse would be present at input 2 of the corresponding trigger. This trigger, 
in switching, would unblock a definite tube in the conductor assembly and, in the circuit of the indicating de- 
vice, a current would be established corresponding to the received code pulse. The circuit would function simi- 
larly to the end of the cycle, at which time the indicating device would have a current corresponding to the re- 
ceived code. If there were no variation in the measured quantity, then no switching would occur in the circuit 
and the magnitude of the current in the indicator would not change. For a new combination of code pulses, a se- 
ties of triggers would, in general, be switched; triggers not taking part in this decoding would be restored to their 
original state by pulses formed at winding W, 5. Here, pulse suppression would not occur, since, in this case, the 
elements of K are not switched. No pulse will appear at input 2 since there will be no pulse from winding W, 4 
of element K to enable sensing in the coincidence circuit. It should be explained that if, in the receiving process, 
the code improperly impinges on the distributor, due either to a transient disturbance in the circuit or to an in- 
acceptably large magnitude of instantaneous noise, the indicating device will not be affected, since the cycle 
repeats itself after 0.2 sec and the actual value of the measured quantity will be restored. 


The functioning of the receiver is shown in Fig. 11. 
The overall error of the design is comprised of the discreteness error, 6, = #{2"-1)] “* (where n is the 
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number of elements in the code) and the summing error, §,, affected by the variation in voltage drop across the 
conductor circuit resistors, which in turn depends on the combination of branches switched. Due to this, the ratio 
of the k'th current i,k to the current i¢ will deviate from the law of geometric progression. The ratio i,/ i,o 
is supposed to equal 2K, i.e., i, = 2, i,%/i,o= 4, .... in-1/i,0= The deviation of magnitude from the 


_ ok 
values defined by this law, relative to the quantity 2" . 1, will give rise to an error 5k) = er st 2 Jae -1), 
2 
where k is the ordinal number of the term in the binary series. 
The currents, i,* and i,0, are expressed by elementary formulas and their computation presents no complexity. 
The maximum summing error for the design will occur when all branches of the conductors are switched: 
n 
5, 2 Sk): In the design considered here, this error amounts to 0.55 %. 
The influence of a series of other factors on the accuracy of measurement, e.g., the inaccuracy of calibrat- 
ing the resistors in the summing conductor circuit, the instability of the voltage source for the conductor circuit, 
etc., may be similarly estimated and held to a minimum. We might mention that it is possible to ignore varia- 


tion in the voltage supplying the equalization circuit if its source is the same as the source of the voltage used in 
the rectification of U,. 


In the pulse-code transmitter and receiver herein described, it seems desirable to replace the ribbon material 
cores in the magnetic distributor with ferrites having rectangular hysteresis loops, in view of their improved tech- 
nology, or to construct the devices wholly of crystal triodes and diodes. 


An important task for further perfection of pulse-code telemetering is the reproduction of digital indications, 
which would allow the avoidance of errors in the receiver's indicator. 


SUMMARY 
1. Pulse-code systems of telemetering possess significant stability under noise, but involve relatively com- acc 
plex apparatus. Therefore, their most felicitous application should be in systems with long, complex links, where tive 
the use of other systems is not feasible or possible, and also in complex automation and remote control units in Wy. 
conjunction with computing devices. ing 
2. The construction of digital reading and registering indicators is the very first problem to be solved in the 
further perfection of pulse-code telemetering systems. 
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RELAY PHENOMENA IN LOOP CIRCUITS CONTAINING MAGNETIC CORES 
WITH RECTANGULAR HYSTERESIS LOOPS 


V. A. Zhozhikashvili and K. G. Mitiushkin 
(Moscow) 


Relay phenomena, observed in circuits containing magnetic cores with 
rectangular hysteresis loops, and containing positive feedback paths, are described. 
An analysis is given of the static characteristics of such circuits, the character of 
their transient processes, caused either by isolated or by repeated disturbances, and 
of the influence of noise pulses. . 


A large group of circuits with magnetic cores is comprised of link circuits, in which a series of identical 
cores are so connected that each core, when switched by an external switching pulse from state 1 to state 0 (Fig. 
1), switches the core following it from state 0 to state 1. In the majority of such circuits, the cores are connected 
according to one or another of the schemes shown in Fig. 2, a and b. Windings Wy, Wz, and Wp are called, respec- 
tively, output, input and supply (shift) windings. A circuit joining two neighboring cores and consisting of windings 
W, and W,, resistor r, and valves B, and B, (Fig. 2, a) or valve By and capacitor C (Fig. 2, b) is called a connect~- 
ing link, or a transfer link. 


Let us assume that an external switching pulse, flowing through 
z - the supply winding, engenders negative ampere turns in the core. The 


action resulting from these ampere turns we will call voiding of the 

core. Voiding will be termed complete if it proceeds along the de- 

-# Ava scending branch of the hysteresis cycle and, during the time of action 
: . of the switching pulse, the core is switched from state 1 to state 0. 


In contradistinction to an external switching pulse, a switching 
pulse entering from the connecting link to the input core winding 


0|-8, engenders in it positive ampere turns, The switching thereby induced 
we call priming of the core. We shall say that the priming is com- 
Fig. 1. Rectangular hysteresis plete if, during the time of action of the switching pulse, the core 
loop. state changes from point 0 to point 1 along the ascending arm of the 


complete hysteresis loop, 


Voiding of a core is accompanied by the appearance of a pulse in its output winding, the power of which 
is the larger, all other conditions being equal, the larger the power in the supply winding pulse, For a sufficiently 
large power in the supply winding pulse, the power in the output winding pulse, during voiding of the core, may 
be significantly greater than the power necessary to prime the following core. In such a regimen, each core is 
a pulse amplifier. If the losses in the connecting link are not too large, then the complete voiding of each core 
is accompanied by the complete priming of the following core, Thus, the process of voiding and priming is trans- 
ferred from core to core, and bears this sustained character for any number of cores inthe chain. This process is 
terminated with the voiding of the final core in the chain. 


If a connecting link is furnished between the last and the first cores comprising the circuit, we shall then 


have a series of connected pulse amplifiers with a positive feedback connection. In the following, we shall call 
such a circuit a looped circuit. 
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A looped circuit with magnetic cores may, if the supply windings are fed by a source of switching pulses, 
possess many stable states. A transition of the circuit from one stable state to another occurs at definite critical 
values of voltages, currents and circuit parameters, and has a snowballing character. 


W, 
($% ¢ 


A 


Fig. 3. A two-half-period looped circuit. 


_ We consider the looped circuit depicted in Fig. 3. This is a development of the circuit shown in Fig. 2, a, 
and henceforth will be referred to as two-half-period. In this circuit, the supply windings of the odd-numbered 


cores are connected in one series, and the supply windings of the even-numbered cores are connected in another 
series. 


We assume that initially the supply voltage U= 0 and that all cores are in state 0, and we call this the 
first stable state. If U~ is now slowly increased, then at the terminals of all the core windings will be observed 
a gradual increase in the amplitude of noise pulses, stemming from non-ideal saturation of the material in the 
neighborhood of point 0, and also from the presence of reverse conductance in the supply circuit valves, Bs and 
B4 which leads to partial priming of the cores. When voltage U,, attains the critical value Uy (Fig. 4) which 
may be called the excitation voltage, the amplitude of the noise pulses becomes sufficiently large to trigger the 
snowballing transfer of the circuit from the first stable state to the second. In this state each source pulse voids 
one of the cores, while complete voiding of a core entails complete priming of the following core. Thus, in the 
second state, each of the cores voids once per cycle, a cycle consisting of a number of half-periods in the supply 
circuit equal to the number of cores in the chain. A transition of a circuit from the first to the second stable 

state without special priming pulses will be called a self-excited circuit. In self-excited circuits, generally speak- 
ing, the amplitude of the noise pulses is attenuated, since a significant portion of the supply voltage is applied in 
voiding the cores. This allows the supply voltage to be significantly increased without disrupting the functioning 
of the circuit. When the supply voltage is made equal to Uyo the amplitude of the noise pulses again attains the 
earlier value and, as a result, a second snowballing process occurs, implying a transition of the circuit to the third 
stable state, characterized by the fact that with each pulse in the supply circuit two cores are simultaneously voided. 
With additional increases in voltage, further skips may be obtained. 


In contradistinction to the transition to the second stable state, the transition of the circuit to the third, fourtl 


- and succeeding states is brought about not only by noise pulses but also by a reverse influence of a voiding core om 
the previous core. 


After the circuit has been transferred to the second stable state by voltage Uy , it can remain in that state 
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even if the supply voltage is significantly reduced. Fig. 4 shows that only at voltage Ue, Will the circuit return 
to its original, nonexcited state. There is thus the typical loop characteristic, peculiar to relay devices. 


The greatest practical interest is presented by the 
Working characteristics of regions I and Il. If the supply voltage 
pulses is chosen to lie within the limits of region I(Uye < U.K 
< Uy), then a self-excited looped circuit is impossible. 
In order that the circuit be transferred to an excited state, 
pulses from an external source are necessary to prime 
one of the cores. Priming may be complete or incom- 
plete. In the latter case, however, the variation in flux 
during priming must be larger than a certain minimum 
quantity necessary for exciting the circuit. As will be 
made clear in the sequel, the length of the transient pro- 
Fig. 4, Looped circuit characteristics. cesses will depend on the magnitude of the flux variation 
OA is the operating pulse curve; OB is the ee 
noise pulse curve; Uy, and Uyo are supply If U., = Uny, then the amplitude of the noise pulse 
circuit voltages, corresponding to transitions for any core is equal to en. After priming one core, the 
to the first and second excited states; U,.) variable signal amplitude for core voiding equals ec. 


is the relaxing voltage, corresponding to If the supply voltage is chosen within the limits 

a circuit transition from an excited to a 

non-excited state; ace the of region Il <U~< Uyg), then the introduction of 

canting supply voltage will produce a self-excited circuit. This 
ay Pe will automatically establish a regimen in which one of 


poise amplitudes ta the opesatng See the cores is found in state 1, and this state will be trans- 
Iand Il; e,3 and@p, are the currents 


out 


corresponding to the noise amplitudes in ferred from core to core, All other cores will remain in 


state 0. If, for example, Uw = Un,, then the available 
signal amplitude will equal e,,, and the amplitude of the 
noise pulses will equal eno. 


Both modes of functioning of the looped circuit, i.e., with self-excitation (region I) and without self -excit- 
ation (region II), may have interesting practical applications. 


Curves OB and OA, bounding the characteristic 
loop respectively below and above, may be obtained in 
the following way. Let one core be wound according to 
the data calculated for the cores in the looped circuit, 
and let its windings, Wy and W2, be loaded as they would 
have to be loaded in the circuit, and then include wind- 

ing Wp together with the valves after 
which incre the supply voltage for e core 
Fig. 5. Different states of a looped circuit, will give 
but with shorted supply circuit valves, then curve OA will 
be obtained. A somewhat higher free-running characteristic may be obtained if the circuit for windings W, and 
W; is open. 


the same cases; U, and U, are values of 


supply voltage. 


Families of loop characteristics, analogous to those shown in Fig. 4, and derived for different circuit para - 
meters, core windings, etc., give a sufficiently clear picture of the static properties and functioning of looped 
circuits. These characteristics are a reflection of the excitation voltage, the releasing voltage, the source pulse 
amplitudes and the noise pulse amplitudes. 


Different stable states may be obtained, not only by varying the magnitude of the supply voltage, but also 
by different values of the circuit parameters. These states differ one from another by the number of voided cores 
during the time of one cycle. Thus, fora very large resistor rt» in the connecting link, none of the cores is primed, 
and the circuit will be in the non-excited state (Fig. 5). With the lowering of resistance tp to a value Tf, one of 
the cores will be primed, as a result of amplification of noise pulses, The circuit will transfer to the first excited 
state, in which each of the cores will be voided once per cycle. This regimen is maintained until t» is lowered 


| 


to a value fo. At this latter value, the circuit transfers to the second excited state, in which each of the cores is 
voided twice per cycle. With further decrease of tp, one may obtain the third, fourth and all succeeding excited 
states, upto the final state, in which each core is voided as many times as the period of source frequency can be 
packed into one working cycle of the circuit. In the case shown in Fig. 5, the circuit contained 10 cores. The 
existence of all stable states has been observed in a circuit containing 30 cores. : 


All states of a circuit may also be obtained by increasing resistor ry from zero. The values of resistor rp at 
which the circuit transfers from one excited state to another when ry is increased do not coincide with the values 
of ty when it is decreased. 


In the case of cores with little noise, self-excitation may, in general, not occur within the attainable limits 
of supply voltage and current, In such cases, self-exitation may be obtained by shunting valves B, and B, with a 
large resistance, thereby causing a partial priming of a core by a pulse of opposite polarity to the supply pulses. 


Special cases of looped circuits with positive feed- 


6, / 2 8 4, back are those circuits containing one or two cores. Such 
ont} circuits may be utilized as magnetic dynamic triggers or 
ng ne% —Tt " relays. Two such circuits are shown in Fig. 6, a and b. 

> he These circuits possess two stable states, and the loop char- 

acteristics, analogous to those given in Fig. 4, may be 

8; +p t Wy w, 3 obtained with them. In the circuit with one core, the 

Z, % : voiding of the core leads to the charging of the capacitor, 
—c> 4; af P Priming of the core occurs as a result of the condenser 
Rp F age discharging through input winding W,. In the circuit with 
“aon, two cores, priming of one core occurs as a consequence 
a of voiding of the other. 


ee ee To function as a control relay, the core is furnished 
with a starting winding W, and an inhibiting winding Wj. The approximate characteristics of a control relay 
for two cores are given in Fig. 7. In comparison with other circuits, the magnetic non-contact relay circuit given 
here possesses both advantages and disadvantages, but an investigation of these questions is not within the scope of 
the present paper. 
An investigation of the transient processes in a 
Tour looped circuit presents significant difficulty in the ge- 


neral case, due to the fact that these circuits contain non- 
linear elements, namely, cores and semiconductors. 


In what follows, certain questions of transient re- 
sponse will be investigated where the rectangular hysteresis 


- Fon loop will initially be assumed to be infinitely narrow 
Frei a Groid (cores with this characteristic will, in the sequel, be called 
ideal), and the direct and reverse resistance of the valves 
Fig. 7. Characteristics of a control will be considered, respectively, as being constant and 
relay. infinite. 


U and I,,,,, are the voltage and current 
in the control path; Uyoigand are 
the voltages in the control path corres- 
ponding to the voiding and releasing re- 


In view of the fact that in the steady state each 
core of the looped circuit must be switched around the 
complete hysteresis loop in the course of one half -period 
of supplying voltage, both for voiding and for priming, 

and investigation of the transient processes in a looped 
, circuit is most felicitously carried out separately for each 
supply half-period. Such an approach was recently employed in investigating the workings of magnetic amplifiers 
(1). 


We shall consider looped circuits with arbitrary numbers of cores. For definiteness , we will take them in 
the form of two-half-period circuits. 


At the input winding W, of some core, in one of the half-periods when there is no current in the supply wind 
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ing of this core we supply a voltage pulse, e, = E, sin wt, of duration t, where t- may have any value from t,=0 
. + oe (T being the period of the supply voltage). This pulse gives rise to a variation A@, of the core 
ux e to 


priming this core for voiding. In the general case, priming may be either incomplete or complete and A® = 26, 
where ®, is the value of the residual magnetic flux in the core, 


In the succeeding half-period, when current begins to flow in the supply winding Wp of the core under con- 
sideration, voiding of the core occurs with a change of flux 4@, during voiding which is generally less than A@_. 
in what follows we consider only the case where A, = A@,,. This will always be the case if the amplitude of the 
emf in the source winding, but neglecting the voltage drop across resistor Rp, and the amplitude of the source 
voltage satisfy the condition 


> 20,. (1) 


The time t, necessary for voiding the core is determined from the following equalities: 


= (1 — = (2) 


During the time of voiding, a pulse is induced in the output winding W, with emf 
Coat = Sin ot = Ew sin ot, 
D 


and length 


For an ideal core, this emf is equal to the voltage approached in the winding W, of the following core, 
inducing a variation of flux in that latter core. This flux variation. 4g, equals: 


— cour = (1 — 208 at). 
0 


Comparing A@, with 4@,, we find 
A®, = Ao = ky. (3) 


Thus, kw shows itself to be the transfer coefficient (amplification factor) of the flux and, depending on 
whether this coefficient is greater than, or less than, unity, there will be, respectively, an increase or a decrease 
in the magnitude of flux variation and in the voiding time for each successive half-period. For ky > 1, the pro- 
cess grows up to the point where the voiding of the core in some half-period produces full priming of the follow - 
ing core (Fig. 8). 


With this latter half-period begins the steady state, during each half-period of which each core in turn is 
voided and the core following is completely primed. For kw < 1, the process, conversely, decays until the flux 
variation during priming becomes equal to zero. 


If we number the cores and the supply half-periods, counting as number one that half~period in which the 
signal was applied to the first core, then the flux variaton in the second core during its second-half-period prim- 
ing will equal A®,= ky 4@4; the flux variation in the third core during the third half-period will equal 


A®, = \ codt = (1 — 008 ate) 
61 g 


Ads = = ki) etc. In general, the flux variation will equal = ky, where n = 1, The core 
voiding time can be determined in the following fashion: 


i 4 
AQ, = (1 — 608 etn) = kw AQ®,, (4) 
from which 
Wy 
costal, = (5) 

In order to determine the number of half-periods 
in the transient process, and the number of the core with 
which complete voiding begins (steady-state), it is neces- 
sary to make the substitution 4, = 26) in Equation 
(4). We then obtain 

kw A A®, 
from which 
1 20 
°8 KO, +1. (6) 


Fig. 8. Process of exciting a looped circuit 
when kw > 1. AB, is the flux variation in 
the first core, induced by signal pulses; 4B, 
and AB; are the flux variations in the first 
and third cores for priming, AB, and AB, the 
same quantities for the second and fourth 


Thus, as follows from (3), the phenomenon of self- 
t excitation is possible only for kw > 1, which is the neces- 
sary condition for the existence of relay-like jumps and stable functioning in several different states of the circuit 
we have been considering. It must hold in all those cases when cores are utilized in commutator, counter, ac- 
cumulator and relay circuits. 


For kw < 1, a self-exciting looped circuit is impossible but, due to the positive feedback, it is possible to 
evoke a significant amplification of noise pulses. In fact, if we assume that in the circuit of Fig. 6, a, during the 
half-period when current is flowing in winding Wp of core II but current is blocked from winding Wp of core I 
by valve Bs, a noise pulse might act in winding W, 4 pulse given by ey = E, sinwt, where t may take any value 
in the range 0< t< T/2. The first such pulse would occasion a partial priming of core I, varying the flux in the 
core by the quantity A@,. In the following half-period, core I would be voided and core II would be primed, the 
flux variation in the latter core being k,,4@p, in accordance with the above. The voiding of core II leads in 
turn to a flux eae magnitude ky 4p. Adding in the action of the noise, we obtain a total flux 
variation in core I of 4g, + kw @n- In the succeeding voiding of core I, a flux change of ky (A @ + kw Sn) 
is induced, etc. 


Continuing this line of reasoning, we do not find it difficult to establish that, in the n'th half-period after 
the appearance of noise pulse en, the flux change induced in core II will be 


Correspondingly, the voiding time in the steady 
state is obtained from Expression (5): 


20, Wpe 
E 


~ 


= A®, (kw + Ky + Rly + Ry), (8) 


where n is an arbitrary odd number. 
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If kw < 1 then, for n = a, we obtain the steady-state value of the flux in core Il: 


kw 


AP hiss = A®, (9) 
The flux change in core II at the (n+ 1) th half-period 
will be equal to 
48; }--------- a A®, = (1 + + 
+ ky + ky +... + 
44°48, a : The steady-state flux in core I equals 
48,... 5, re 
Thus, the steady-state flux a the resulting 
Fig. 9. The process of trigger ex - flux change is proportional to the flux induced by the noise 
citation under the influence of noise pulse but, in the first core, is larger than this initial flux 
pulses for ky, < 1. By and By, are, 1 
respectively, the inductions in cores (1- ke)’ w) 
I and Il, and AB, are the variations the noise pulses taken off the output winding of core I, as 
of induction with each noise pulse. compared to the same pulses in the absence of positive 


feedback. Fig. 9 shows the variation of magnetic induc- 
tion in cores I and Il for kw < 1, corresponding to the process just described. 


In circuits consisting of real magnetic cores, having a coercive force Hj 0, the transition from the 
non-excited state to the excited state may occur only if the signal pulse e, or the noise pulse e€, exceeds a de- 
finite threshold value. 


In order to determine the aforementioned threshold value and to elucidate the conditions under which it is 
possible for the circuit to function in excited states with completely voided cores, we shall consider not infinitely 
narrow hysteresis loops, but rectangular loops with a finite width 2Hj. Here we will take the circuit to be of the 
form given in Fig. 3. 

we. ws 

We introduce the following notation; a= te tw , b= S-10™, where S is the cross-sectional area of 
» where L is the average length of magnetic line, 
D,= + hy), + kw) + Ds, 


the 
core, 


we 


AB, 
where AB, and AB, are the total variation in inductance during the voiding and the priming of the core, respectively. 
For a sinusoidal source of pulses, the process of switching cores may be thus described: 
= Asin ot — — 


The process of varying the induction during priming, under the same conditions, may be written as: 


ue kw Asin wt — Dy. 


—coswt,) — Dit, = AB,, 


A 
(1 —coswt,) — = (13) 


where t, is the time during which the induction varies by the amount AB, during voiding of the core. 
As in the case of the ideal cores, we compare Expressions (12) and (13), obtaining 


AB, = kw AB, — D,t, 


kp = kw — Ds (14) 


Thus, in contradistinction to the circuit with ideal cores, in the circuit with real cores, kp y kw. The value 
of kg depends on the ratio of the switching time to the induction variation during switching. The condition ky > 1 
{s no longer sufficient to guarantee the possibility of jumps, and of the circuit functioning in excited states. 


Reasoning the same way as for the case of ideal cores, it is possible to state that if in one of the half-periods 
when current is absent from the source winding of the first core, priming of this core is occasioned by a pulse e, 
of magnitude 4B., then in the following half-period, when the first core voids, the flux variation AB, will in 
general be less than or equal to. AB... As in the ideal case, we will consider only the situation when AB, = ABe, 


For this, the following condition must hold: ‘ 
a 
2A T 
= — Diz > 2Bo. 
We assume that the time of voiding t, is determined by Expression (12). Let ty = fy(AB,). Voiding of wr 
the first core entails priming of the second core, wherein the induction variation is determined from (13): wat 
th 

ABy = ky AB, — = kwAB, — Dif (AB,). 
In the second core, during the following half-period, voiding time is determined from the equation e 
A (4 — Dyty = AB,, 

sno 

pes 

ty = (ABs) = [kw AB, — (A4B,)) = (AB,). 
When the second core voids, the induction variation in the third core equals ‘an 
AB, = ky AB, — = AB, — ky Dy f, (AB;) — Daf (AB,). 
Reasoning similarly, we may write for the voiding of the n*th core during the n‘th half-period: 
tha 
AB, = [1 — eof, (AB;)] — (AB}). As, 


The induction variation in the (n+ 1)th core during the n'th half-period may be determined from the 


pression: 
ABais = AB, — Daf, (AB,) — Daf, (AB,) —.. 
hy Dsfn—1 — (AP,). 


| Integrating, we get 
on 

| | 
| or 

4 
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Thus, although the analytic investigation of the transient processes is theoretically possible, in practice, it 
involves considerable difficulty. 


NIN & 


48, 26, 


48.,i,48, 48, 24, Fig. 11. Excitation process of a 
looped circuit for sinusoidal switch- 
Fig. 10. Consturction of the curve, ing pulses, and Hj 0. 
= f(AB). 


A graphic determination of the conditions for exciting the circuit, and of the length of the transient processes, 
may be carried out as follows. From Equation (12), the dependence of t on AB is constructed, as shown in Fig. 10; 
above this curve, using it and Equation (14), one constructs the dependence of kp and 4B, Functioning in an ex- 
cited state is possible only on the condition that, for values of AB lying in the interval B,;, < 4B < 2Bp, kp assumes 
values greater than unity. The value of B,,j, for which kp = 1 thus defines the minimum signal necessary for ex- 
citation of the circuit. If the magnitude of the signal is known, then the length of the transient process may be 
determined from the curve of kg. We do this by determining AB, from the known signal amplitude. Using this 
value of AB,, we obtain kg, from the curve, Muftiplying AB, by kg,, we obtain a new value AB, and so on, until 
the point is reached at which AB = 2Bo. 

In Fig. 11, the region of relay action is shaded. On the axes of the abscissas and the ordinates lie the induc~- | 
tion values for, respectively, even and odd numbers of cores. 


SUMMARY 


1. Looped circuits with magnetic cores possess, generally speaking, many stable states wherein the transi - 
tion of the circuit from one stable state to another occurs for definite values of current and voltage, and has a 
snowballing character. This phenomena stems from the amplifying action of the cores and from the presence of 
positive feedback. 


2. The stable states, corresponding to regions I and II of Fig. 4, may be used to obtain counting -switching 
devices, working under a regimen either of external excitation (region I) or of self-excitation (region Il). 


3. A looped circuit consisting of ideal cores (Hj = 0), with the condition that kw > 1, may not, in principle, 


remain in a non-excited state, since an arbitrarily small disturbance suffices to transfer the circuit to an excited 
state. 


4. A looped circuit consisting of real cores (H; = 0) may be taken out of the non-excited state providing 
that the disturbing action of signal or noise produces a switching of a core which exceeds a certain threshold value, 
AB min: 
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DYNAMIC CHARACTERISTICS OF CORES WITH RECTANGULAR STATIC 


HYSTERESIS LOOPS (INFLUENCE OF EDDY CURRENTS) 


M. A. Rozenblat 
(Moscow) 


This paper deals with the influence of eddy currents on the form of the dynamic 
hysteresis loop, the magnitude of the differential magnetic permeability and the magnitude 
of the dynamic coercive force of cores with rectangular static hysteresis loops. Analytic re- 
presentations are found for the dynamic hysteresis loop for sinusoidally varying induction, 
for sinusoidally varying field intensity, and for dc magnetization of the core. 


Experimental data are cited in confirmation of the theoretical results. 


INTRODUCTION 


For the cores of magnetic amplifiers, widespread use is made of alloys with rectangular hysteresis loops (e.g. 
5ONP, 65NP, N34K29MZP, etc.), which makes possible significant improvement in linearity, increase in gain, and 
decrease in lag of amplifiers with positive feedback [1]. Such alloys have also found wide application for other 
electromagnetic devices, for example, in magnetostatic triggers, in memory elements, and in pulse transformers, 


In the function of such devices, the width and form of the dynamic hysteresis loop are significant. For ex- 
ample, the magnitude of the differential magnetic permeability, fp = dB/ dH, on the “vertical” portions of the 
dynamic hysteresis loop determines the value of the gain for a number of types of magnetic amplifiers with posi- 


tive feedback. 


a 


4, 


— 
W(-8,,0) 


Fig. 1 
N34K29MZP, the ratio B,/B, attains a value of 0.95-0.99, and up a value of 10’ gauss/oersted. 


In a dynamic regime, the cores and their lamina or turns are magnetized inhomogeneously across their cross 
section. Therefore, we shall understand by the dynamic characteristics of a core the average dependence, across 


In order to appraise the quality of cores with magnetic hysteresis loops, 
one frequently checks their dynamic characteristics experimentally (e.g., 
by means of an oscillograph). To interpret the measurement results, it is 
necessary to know the dependence of the dynamic characteristics on the con- 
ditions of measurement. 


The aim of this paper is the theoretical investigation of the influence 
of eddy currents on the dynamic characteristics of toroidal cores with rect- 
angular static hysteresis loops under different conditions of magnetization. 
We ignore the influence of magnetic viscosity, For strip material about 0.03 
to 0.04 mm and above in thickness, this influence is ordinarily small in com- 
parison with the effect of eddy currents [2]. It is assumed that in the limit 
the static hysteresis loop has the ideal rectangular shape (Fig. 1), characterized 
by the residual inductance B, being equal to the saturation induction B, 
and the differential permeability 4 being equal to infinity for values of 
| < 
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its cross section, of the magnetic induction on the strength of the field induced by the current in the magnetizing 
winding. In the present work, we assume that the ratio of the inner and outer core diameters is close to unity, and 
we therefore ignore the inhomogeneity of the core magnetization, due to the lengthening of the magnetic path, 
and the corresponding decrease in field intensity, stemming from the increase of diameter in the portion of the 
core under cofisideration. 


Equation of the Dynamic Hysteresis Loop. 
We consider the circuit of Fig. 2, for which we have 


w Ri =e, (1) 


where e is the voltage necessary to switch the core. 


The influence of the eddy currents is studied by 
means of a method due to V. K. Arkad‘ev [3, 4]. For 
this purpose, we consider the core to be in the form of 
a lamina with length L, height b and width d (Fig. 3). 
For the strips (windings) of the core, L = nj ,where n is 
the number of tums about the core, is the length of 
one turn, b is the core depth and d is the width of the 
strip; asa rule, L» b» d. 


We suppose that the initial magnetic state of the 
core corresponds to the point N (—B,, 0) of the hysteresis 
loop (Fig. 1). If the field intensity, induced in the core 
winding w by the current i, exceeds the coercive force 
b Hg, then the induction in an infinitesimal strip of width 
dz is changed instantaneously from ~B, to +B,, The flux 
change in the lamina begins at the circumference and 
proceeds inward. 


Let the flux penetrate to depth z. Then, further 
variation of the induction in the two strips of width z 
zilaz \iz (Fig. 3) induces an emf (taking d« b), e = —2B,° 
de yy -2 10%, tna strip of width 2, this emf induces an 


‘ eddy current, if = eY, where Y is the conductivity of this 
strip, the magnitude of which is found from the formula 


= 424. , being the specific conductivity. 


Fig. 3 This current induces a field intensity 


= —0.8nB, 72 = 10", (2) 


During the entire process of switching the core, for the center of the laminar cross section (z = d/ 2), in which 
the induction still has not varied, we have 


H+ He = H, (3) 
where 


= (4) 


By taking into account the fact that the core flux @ is related to the quantity z by the formula ¢ 


we obtain from Equations (1-4), 


M (® + ®,)d® + Nd® = (h—H,) dt, (5) 
where d 
0.8nB, yd? 
0.Arw* (7) 
N = 
and 
h= (8) 
The core begins to be switched only after h, in changing from zero, attains the value H, at time t= t' 
{h(t") = H,]. Thus, integrating Equation (5) from t’ to t, which corresponds to the flux change from—@, to @, . 
we obtain > 
t 
or, using the fact that at all times @+ @, = 0: 
N /m 2 
Substituting Expression (10) for the flux in (1), and solving the resulting equation for the current, we find 
w (h — H,) 
=R 
RM.-108 
For the field intensity (4), we have 
N 
(h— H,) 
H= h— — 
Ne 
From (10) and (11) we obtain the equation for the ascending arms of the dynamic hysteresis loop: 
M 
+P h an 
@+0,) 


| 


H= (12°) 


where § is the core cross sectional area, and B = @/S is the average or apparent induction. 


As might be expected for y = 0 (which occures with ferrites), and also for d ->0, we get H = H,, and the 
dynamic and static hysteresis loops coincide. 


If, in Formulas (12) and (12"), we set, respectively, @ and B equal to zero, we then obtain the value of the 
coercive force for the dynamic hysteresis loop: 


H, +79, h 
(13) 
1 +70, 


For h = const. = Ho, we see from (13) that H,p = Hg. 
We now consider certain examples in which the formulas so far obtained are applied. 


Sinusoidally Varying Induction 


Let e = E,,sinwt. The induction will vary sinusoidally for R= 0. Integrating Equation (1) for this case, 
and bearing in mind that at t = 0 the indiction B = —B,, we obtain 


En 
B= —B,+ ows 10% (1 — cos wt). (14) 
We set 
En 


h Em 
For R=0,7-=0, and = "10° sina. 


Therefore we have from (12) , using (15): 
H = H.+ sin at. (16) 
From (14) and (15) it is not difficult to obtain 


sinwt = (+ 


and Equation (16) for the dynamic hysteresis loop takes the following form: 


Setting B = 0 in (17), we obtain the following representation for the dynamic coercive forces 


or 
69 
% 


p = 3H,. The dotted lines indicate the static loop. It is obvious from Formula (17) and Fig. 4 that the field 


Fig. 4 shows the construction of the dynamic hysteresis loop from Formula (17) for the values k = 1 and 


intensity H is a nonlinear function of B during the core switching 


8 process and, consequently, that the differential permeability is a 
variable quantity. From (17) the following expression for the different- 
: : ial permeability is easily obtained: 
2k (4 fi 1 BY 
dH emo? By) _ BY 
M se (1+ 
For B = et. k— 1)B, , Up = ©, and the field intensity H. 
attains its maximum value, equal to 
nu 
Hm =H. + or, fe 
an 
i? For further increase in B, the field intensity H is decreased, 
and fp assumes negative values. When k = 1, corresponding to the induction varying from —B, to +B, during a 
half-period of the supply voltage, H attains its maximum, and fp = © for B = +0.5B, (Fig. 4). If k > 4/3, then 
H does not attain its maximum, since this would have been observed at B > Bs. 
be 
inusoidally Var ield Intensi 
We let for 
0.4nwE,,, 
h= sin wt = H»sinat. 
We obtain the same law of variation for H if R-®on, but the ratio E,,/R= const. In this case N->0 and Dc 
Equation (10) assumes the following form: 
mil 
© +0, = [Hm (008 wt’ — cos wt) — (t—t’)], (19) For 
and, from (12), we obtain by | 
whe 
H=h=H,sinot. (20) con 


Since H,,sinwt'= Ho, we can get rid or t and t’ in Fomula (19), obtaining for the dynamic hysteresis loop 


— H, (arcsin — are sin 


(21) 


The case when Hy, » H, is particularly interesting. 
Since, in this case, the switching of the core occurs for a value H « Hy), we may substitute in Formula (2!) 


| 

| 
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Making these substitutions in (21), we find * 


H—H 
B = By 
Thus, for Hm » H, and also for a linear variation of H in time, the induction B varies quite proportionately 


to H, but the differential permeability remains essentially constant during the entire process of core switching, and 
equals 


It might be mentioned that fp is frequently measured by oscillographic methods in order to determine the 
number of cores to be used in magnetic amplifiers [1]. For this purpose, the cores are magnetized by the sinusoidal 
field (20), for which Hj, » Hc. Then, the sine occuring in Formula (20) may be replaced by its argument (20°), 
and for the emf induced i.. core winding w, we obtain 


dB dif 
e= — wS10* — w510 = — owSH 
It is clear from Formula (23) that the variation of 
be accounted for in carrying out the measurements. 


If we set B = 0 in Formula (22), we obtain the following expression for the magnitude of the dynamic coercive 
force when Hp, » Hep: 


He +O, VOM Am = He + VIO (24) 


DC Magnetization 


Let e = E = const and, correspondingly, h = Hy = const. In this case the dynamic hysteresis loop is deter- 
mined from Formula (12) or from Formula (12"), and the dynamic coercive force, after substituting Hy for h, from 
Formula (13). It is easily seen that for Hy = H,, we have Hp-p = Hg. 


In many types of magnetic amplifiers with positive feedback, the magnitude of the load voltage is determined 
by the value of the induction which is established in the core during a half-period'(T/ 2) of the supply voltage, 
when acted upon by a de signal (cf., for example, [1], par. 7-5). Taking into account that t'= 0 in the case under 
consideration, and substituting t = T/2= 1/ 2f, we find from Formula (10) 


N i 
By = — V + Ho— (25) 


The gain, either of voltage or of current, is proportional to the derivative [1] : 


28 Y (He — (26) 


ei does not depend on the magnitude of 


* It might be mentioned that the condition H « H,, corresponds to a linear variation in time of the field intensity 


H=h=H,et (20") 
and Formula (22) may be obtained directly from Equations (10) and (20). 
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the signal, Hg (or E). The presence of eddy currents not only decreases the gain, but also makes it a function of 
the signal. An increase of the latter causes the derivative, and hence the gain, to decrease. 


If in Formula (25) we set By = Bs, we obtain the following expression for Hp maxs the field intensity of the 
signal necessary to magnetize a core completely in one half-period of supply voltage: 


Hymax= He + 4/MS*B,? + 4/NSB,. 


Thus, the presence of eddy currents entails the necessity of increasing the field intensity of the signal by x, 
the quantity 
AH, = = 
It is also possible to express AHy by means of the dynamic coercive force under sinusoidal induction. Indeed, 
taking into account Equation (18), we find that for k = 1 
cu 
AH, = 0,2 (fl .p— H-). 
Thus, if Hg is known, it is sufficient to measure Hep under sinusoidal induction for evaluating the influence cE 
of eddy currents on the behavior of a magnetic amplifier. an 
the 
Influence of Short-Circuited Windings and Measuring Circuits on the Experimental 
Determination of Dynamic Characteristics 
In certain cases, the presence of shorted core windings, and also the presence of finite resistance in the measur- 
ing device, can significantly distort the results in measuring the parameters of a dynamic hysteresis loop. 
We suppose that around the core (in addition to the basic magnetizing winding, w) is a second winding w,, 
sharted by a resistance Ry. Then, on the left side of Equation (3), it is necessary to append the field intensity 
Hk induced by the current i, in winding wy. For H, we have 
0.4nw?, 
k= Ri, 1-108’ 
we obtain, instead of (9), the following expression: 
of 
cha 
line 
and the equation of the apparent dynamic hysteresis loop takes the following form: (22) 
K M 
Wh+ eon 
It may be concluded from (29) that shorted windings have a completely different effect on the form of the pert 


dynamic hysteresis loop than do eddy currents. Therefore, in calculating the influence of the latter, no recourse 
was had to introducing “equivalent” shorted windings, as is sometimes done. 


If the induction is varying sinusoidally (14), it is possible to obtain the following expression for the 


ay 
| 
* Si 
the | 
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apparent dynamic hysteresis loop, using the previously investigated methods: 


H =H, +0®,[K + Mo, (1+ — (1 +#)- (30) 


Letting B = 0, we find from (30) that the presence of a winding w, shorted by a finite resistance, leads to 
an increase in the apparent dynamic coercive force by the amount 


AH, = KY 2k —1. (31) 


EXPERIMENTAL 


In a work by the author [4] it was shown that the method herein employed for studying the influence of eddy 
currents permits of sufficient accuracy in computing the process of magnetizing strip cores of alloy N34K 29MZP 
with rectangular static hysteresis loops. Here, we cite certain results of investigations of the dynamic character- 
istics of toroidal cores made from alloy 65NP. The cores were made up of rings 0.35 mm thick, 30 mm outside 
diameter and 20 mm inside diameter. The specific electrical conductivity of the alloy equals 4 x 10* (ohm cm)“, 
and the static coercive force of the cores was Hp = 0.06 oersted. In evaluation of the formulas, instead of Bs, 

the value of the residual induction By was used; this value turned out to be 12,000 gauss. 


8 8 


Fig. 5 Fig. 6 Fig. 7 


Figures 5 and 6 give oscillograms of the dynamic hysteresis loops, taken off, respectively, for frequencies 
of 50 and 500 cps, with sinusoidally varying induction (k = 1). It is clear that the general, somewhat unusual, 
character of the dependence B = f(H) corresponds to the theoretical curve of Fig. 4. The oscillogram of Fig. 7, 
obtained with a sinusoidally varying field intensity, shows that in this case the dependence B = f(H) is virtually 


linear during the course of practically the entire magnetizing processes of the core, in accordance with Formula 
(22). 


Table 1 contains the results of computations using Formula (18), and the results of measuring the dynamic 
coercive force Hep for various values of k and of frequency f with a sinusoidally varying induction.* Table 2 
gives the results of computation based on Formula (24) and the results of measuring Hep with sinusoidally vary- 
ing field intensity. 


In Table 3 are given the results of evaluating Formula (23) and the results of meauring the differential 
permeability pp for a field intensity varying sinusoidally with a frequency of 50 cps. 


The results given in Table 1-8 show a satisfactory agreement between the computed and the measured values 


* Strictly speaking, induction varies sinusoidally only for k = 1. For > 1, the induction varies from ~B, to +B, 


according to the cosine law, remains invariant during a certain time interval, and then varies again according to 
the cosine law from +Bs toBs. 


TABLE 1 of Hp and #p; however, it should be borne in 
mind that, in measuring these quantities by the 
Hep. Cormed oscillogram method, one obtains an accuracy no 


computed better than 5-10 


The same agreement between theory and ex- 


‘ 50 0.36 0.355 periment was also obtained with cores prepared from 
2 50 0.63 0.57 0.133 mm thick 65NP alloy. For cores of material 
3 0.05 and 0.02 mm thick, satisfactory agreement 

: 50 4.08 0.94 between computed and measured characteristics 

6 50 1,20 : 7 was obtained only with magnetization by sinusoidally 
ii ee en varying, or by constant, fields, with Hm Or Ho» H,. 


For values of Hy very close to H,, and with sinusoi- 
dal induction, the measured values of H,p turned 
out to be significantly higher than the computed 

ones, which might have resulted either from the influence of magnetic viscosity, or from the inhomogeneity of 

magnetic properties along the cross section of the material. 


TABLE 2 TABLE 3 
| Hep» Oersted | uD» gauss/oersted 
Hm, | computed Hy», 0erstedjicomputed 

cmeme measured fom (24) from (23) | measured 
2 50 0.75 0.82 2 15 700 17 000 

4 50 1.07 1.14 4 11 100 11 800 

6 50 1.46 1.38 6 9 050 8 600 

8 50 1.62 1.56 8 7 850 7 700 

10 50 1.82 1.77 10 7000 6 800 

12 50 2.06 1.94 12 6 400 6 000 

SUMMARY 


1, The influence of eddy currents on the processes of magnetizing cores with rectangular static hysteresis 
loops may be studied by a method due to V. K. Arkad'ev. 


2. The form of the dynamic hysteresis loop and the magnitude of the dynamic coercive force depend, in 
an essential fashion, on the law of variation of the magnetic induction or of the field intensity. — 


3. With sinusoidally varying induction.it is possible to observe a peculiar form of the dynamic hysteresis 
loop, characterized by the following behavior; as the induction increases, the field intensity increases to a certain 
maximum value, at which the differential permeability becomes infinite, then decreases again. This allows the 
differential permeability to assume negative values. 


4. Eddy currents cause a lowering of the gain and a disturbance of the linearity characteristics of magnetic 
amplifiers with positive feedback, controlled by dc signals. The influence of eddy currents may be estimated by 
measuring the coercive force with direct current, and sinusoidally varying induction. 


5. The effect of shorted windings is different from that of eddy currents. It is impossible to simulate the 
effect of eddy currents on the dynamic hysteresis loop by the introduction of “equivalent” shorted windings. 


6. A significant influence of eddy currents of dynamic hysteresis loops is disclosed at frequencies as low as 
50 cps, even for core materials 0.05-0.1 mm thick, and for material thinner than that, particularly with a sinusoid- 
ally varying field intensity. 
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SOME OPTIMAL RELATIONSHIPS IN IDEAL AC-CONTROLLED 
MAGNETIC AMPLIFIERS® 


K. S. Volchkov 
(Moscow) 


The functioning of ideal ac-controlled, active-load, choked magnetic amplifiers 
is investigated. Optimal relationships are deduced for amplifiers designed either for 
amplifying signals of one fixed frequency or for amplifying ac signals in some given 
frequency band; such relationships guarantee maximum gain. The connection is esta- 
blished between an amplifier's frequency distortion coefficient and its time constant. 


At the present time there is a widespread use of semiconductor amplifiers without vacuum tubes for the 
amplification of ac signals of audible and higher frequencies. However, for ac amplification, it is sometimes more 
suitable to use not crystal, but magnetic amplifiers, since the latter, in certain respects, possess more favorable 
properties. Thus, for example, magnetic amplifiers have lower excitation thresholds than do crystal ones, they are 
more easily impedance matched with arbitrary low-resistance transducers and receivers, allow more successful amp- 
lification of very low frequencies, and do not have limitations on their output power. Moreover, magnetic ampli- 
fiers show little variation in parameters from one replicate to another, and their characteristics are quite stable 
under variations in temperature. Considering that the properties of crystal and magnetic amplifiers frequently com- 
plement one another, it sometimes appears reasonable to use the two types in conjunction. 


The present paper is devoted to an exposition of the optimal relationships for ac-controlled magnetic amp- 
lifiers, Herein are submitted to analysis ideal choking amplifiers with active loads, whose ac output signals are 
transformed into the modulating vibrations of a carrier frequency. If two half-period rectifiers are used, this fre- 
quency is twice that of the amplifier supply. In the sequel, such amplifiers will be called ac amplifiers. 


Let the transducer signals entering the magnetic amplifier input (Fig. 1) induce an emf, the angular frequency 
of which is Q, with amplitude Ep,,. (For simplicity we shall assume that the signal transducer has only an active 
resistance, Rp, although analogous reasoning could be followed even for a transducer with complex impedance.) 


In the control circuit of the amplifier will flow a current, the instantaneous value of which is determined 
from the following formula [1}: 


iy = Rd sin (Qt — g), (1) 


where & = Rp/ Ry is a coefficient characterizing the ratio between the transducer's active resistance and the con- 
trol winding resistance, and 7 is the time constant of the control circuit. 


Within the limits of the operating region of its characteristics, 1; = f(I,) (where I}; is the load circuit current), 
the current gain of an ideal amplifier, and its time constant + depend neither on the magnitude nor on the fre- 
quency of the signal. Consequently, there will correspond to the control current iy a completely determinate 
load current,®*the variable component of which may be found from the current gain Kj by the use of 


* Pesented at a seminar on magnetic amplifiers at IAT AN SSSR, Feb. 13, 1957. 
** This take the form of the average current value during a half-period of supply voltage. 
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(2) 


and 


tive 


obt. 

whe 
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the following equation: 


pm 


The effective value of the power thus dissipated by the active load resistance R,; will equal 


P, = iP, = Ru. 


Substituting here the effective value ly frora Expression (1), and the known representation for K, [1], we 
obtain 


+ ay TF (2) 


08 Wy and W.~ are the number of turns in the control and ac current windings, and Kp, is the feedback coeffi- 
ent. 


Fig. 1. Circuit connecting amplifier and signal 
transducer, 


(2) in the following form: 


4fn 
@) 


where f is the supply voltage frequency, and 9 is the efficiency of the magnetic amplifier. 
Expressing the time constant of the control circuit. r in terms of the time constant fy: of the amplifier 


itself 
ty 
tee (4) 
and substituting this in (3), we get 
4fn ty 
(5) 
Ry Ry a+ 


We find the optimal value of ry (given Q, Ry, @,etc.), for which P,, attains a maximum. Taking the deriva- 
tive, 3 P},/d Ty» and equating it to zero, we find 


i 


11 


* 


Using (4), this relationship may be put in a simpler form: 


i 
Topt == (7) 


It follows from this equation that increasing the frequency entails choosing a lower value for the time con- 
stant. 


In evaluating and comparing magnetic amplifiers, it is convenient to use the power gain coefficient. How- iz 
ever, for ac magnetic amplifiers the power gain coefficient may be expressed in different ways, and an established 
method of defining it does not exist. Certain authors understand by “power gain coefficient" the ratio of the amplifier 


Output power to the apparent input power: 


(Al,)*Ry 
(Al,R, V 1 + 


>, = 


If the given equation takes the form 
4fn 


then it may be seen that the maximum gain coefficient, Kp,, is reached for Ty —?>o., 


Despite the continous increase of K with increasing ry, the available power in the load, beyond a certain 
optimal value of ry, determined by Equation (6), will decrease. This is explained by the fact that, with increas- 
ing ry, the apparent power at the amplifier input decreases more slowly than the output power. It is for this reason 
that the coefficient K,. has, in our opinion, formal value and practical inconvenience, since it indicates only the 
relationship of amplifier input and output power, and does not characterize the available power at the load. 


We now consider the influence of ry on the magnitude of the apparent power at the amplifier input. From 
the expression for the apparent power. 


tag 
it can be found that, if «+ 1 < ¥2, then the maximum induced power in the control winding occurs at Ty =% and 
if a+1> ¥2, the maximum power is attained with time constant | ita 
for 
i 
ty =a V (1+ a)? —2. (10) 
A graphic representation of the function Sy = f(ry) for 2 = 316 is given in Fig. 2. Obviously, the func- 
tion Sy = f(Q). will have the same form. 
Since it is clear that a + 1> V2, the transducer signals are computed necessarily on the basis of maximum 
amplifier output, obtained by satisfying Condition (10) and determined from the expression 
in 
1 
Comparing Expression (6) and (10), one may see that after the value 7,, determined by Equation (10), fs ‘dis 


attained, the magnitude of the apparent power at the amplifier input will fall, while the power in the load wind- 
ing will continue to increase until such time as ry attains the value determined by Relationship (6). A= 


18 


ir 
= 


The power gain coefficient may also be expressed as the ratio of output power to the active power dissipated 
in the control winding: 


Al, \2 Ry 
But in this case, the ac gain Kp. turns out to equal the de gain Kg and, just as for K,., does not character- 
ize the amplifier output power, since the admission of active power at its input depends on the signal frequency. 
sy 
ds 
| 
a-7 
20 
15} 
a=3 
as 
4 = 
a / J 6 7 
sec 


Fig. 2. Dependence of a magnetic amplifier's apparent in- 
put power on its time constant. 


It is always satisfactory to define the power gain coefficient as the ratio of the effective power produced in 
the load to the maximum power which the signal transducer is capable of providing. Such a definition is advan~ 
tageous in this respect, that, by means of the gain coefficient, one may judge the increase of effective load power, 
and the degree of utilization of all amplifying devices in units comprised also of signal transducers, Moreover, 
it appears to be sufficiently general, since it is acceptable for arbitrary amplifier types, including those amplifiers 
for which the signal transducers are not generators, but require energy (e.g., Ramey magnetic amplifiers [2]). 


The maximum power which a transducer is potentially capable of providing will equal 


Py= (12) 


D 


Dividing the left and right member of Equation (5) by the corresponding members of (12), and replacing R, 
in the latter by Rp/ a, we obtain an expression for the power gain coefficient: 


at, (13) 


= 


where A = = Fig. 3 give graphs of Function (13) as a function of Ty for different values of Q, and with 


A= 10° and a = 3. It is apparent from the curves that K,, has a clearly defined maximum which occurs when 
Condition (6) is satisfied. The greatest possible gain coefficient for the given signal frequency will be called the 


optimal power gain coefficient, Kopt- Substituting the value of Q from (6) in Equation (13), we find 


= (i+-a) ° (4) 


In accordance with Equation (14), all the optimal gain coefficients on Fig. 3 fall on the straight line OA. 
If, in Equation (13), we set Q = 0, we then obtain the expression for the dc gain coefficient; 


at, 4K 9% 


Kn = (25) 


In Fig. 3, this last expression corresponds to the line 2 = 0, the tangent from the origin to the family of 
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Fig. 3. Dependence of power gain coefficient Kp, om the 
magnetic amplifier's time constant, for various values of 
input signal frequency. 
Fig. 4 represents the same Function (13) as it depends on 2, with ry as a parameter, in the case where 
a = 5. The graphs show two families of curves, corresponding to two values of the coefficient A in Equation (13). 


In order to satisfy Condition (6), Equation (13) is transformed to; 


from which it is easy to see that the geometric locus of the optimal value of K, for variable frequency Q isa 
hyperbola. In Fig. 4, these hyperbolas are represented by dotted lines for each family of curves. The curve Kp, = 
= f(Q) is then the hyperbola made up of the points for which Equation (6) is satisfied. am 


The gain coefficient Kpy expressed by Equation (13), will have extremal values for variable a and fixed 
Ty (Fig. 5). To determine the values of a for which Kp, attains a maximum, we set the derivative oKp,/ a 


an) 


| 
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Fig. 4. Frequency characteristics of a magnetic amplifier for various 
values of its time constant. The dotted lines represent Kopt = £(Q). 


150 


Fig. 5. Dependence of the gain coefficient. on the ratio of the impedances 
of the transducer and of the amplifier control (@ = Rp/ Ry) for ry = 10" 


and A = 10° 
Equation (17) essentially expresses the condition that the transducer impedance be matched to the magnetic 
amplifier input impedance zy since, in this case, 
Ry = Ry Vi = (18) 


It is necessary to determine from which of the two mutually exclusive conditions, (6) and (18), the greatest 
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value of the gain coefficient is obtained. In satisfying Condition (18), the value of the gain coefficient will equal 


Kon ke (19) 


a 


Fig. 6. Dependence of the gain coefficient Kp, on the magnetic 
amplifier parameters. W/W. (solid line) and Ry;/ Ry (dotted line). 


The condition introduced above, (7), for obtaining the optimal power gain coefficient allows a rational 
choice of parameters for an ac magnetic amplifier. If we substitute in Equation (7) the value 


t= Wy (21) 
WR, (1 +a) 4fn(1— Key) 


then we can get the expression 
n F W,R, 
R, — Kg) +a) (22) 


establishing the connection between the basic amplifier parameters necessary for obtaining optimal amplification 

for the given signal frequency, F = Q/ 2x. For example, using the charging coefficients Ksy and K3., of the com- 
trol and ac windings, plus the cross sectional areas of their conductors, qy and qu, we can from Equation (22) 
the required ratio of their aperture areas: 


2 


(23) 


R, W. F 


Fig. 6 indicates the influence of the ratios W)/W,, and Ryy/ Ry on the gain coefficient Kp, for a = 3, 
n = 0.8 and Kfp = 0.95. 


gai 


or 
Considering Relationships (17) and (19), we find 
K. 
| It follows from this that Kop, > K ‘op, » and that for increasing @ the difference between them decreases, 
att 
100) 
Lf qu 
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The functional dependences introduced above were obtained as applicable to magnetic amplifiers of power 
or of current. For the use of magnetic amplifiers as voltage amplifiers, it is also possible to find the most advanta- 
geous amplifier parameters. Expressing the voltage gain coefficient in the form 


AI,Ry 
AI,R, V1 + 


and using Relationship (21), we get 


K 


(24) 


By using (24), we find the turns ratio of the control and ac windings for which maximum amplification is 


attained for the given signal frequency: 


All the previous results dealt with the case when the signal frequency remained fixed, For the behavior of 
an amplifier in a certain frequency band, the computations may be carried out starting either from a given fre- 
quency distortion or from the maximum average power gain coefficient.* The average gain coefficient in the 
given frequency band from Q, to Q, (Fig. 4) is found from the following expression: 


Aaty 


Qs 
1 
Krav = Kon (0) 40 = 


Setting the derivative @Kpy gy/@ ry equal to zero, we find the condition for which the average power 
gain coefficient in the given band will be a maximum: 


i+a 


For amplifier computations starting from the frequency distortion, there is usually given a coefficient of 
frequency distortion. Such a coefficient (for example, the voltage at frequency Q, with respect to that at 9,) may 


be found from Equation (1): 
=| V 1 + (28) 
1+ 


(29) 


Analogously, if the power coefficient of iB, enasiten is given, we have 


M,= 


Solving this equation for 7, we get 


t= (30) 
Pp 


* It is clear that in those cases for which the signal always has one and the same non-sinusoidal form, the amplifier 
computations must be carried out with this peculiarity taken into account. 


| 
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Tests carried out on a series of magnetic amplifiers, with Kp, =< 0.95, showed that Equations (29) and (30) 
are in sufficiently good agreement with experimental data. This allows one to assume that the time constants of 
real ac magnetic amplifiers remain practically invariant with changes in signal frequency (of the order of 1/5 
of the supply frequency). The value of Formulas (29)and(30) consists not only in the possibility of computing from 
them the time constant, given the frequency distortion, but also in the possibility of determining r by experimental 
determination of the frequency characteristics. 


It is clear from Fig. 4 that in those cases where it is necessary to have the very lowest frequency distortion, 
it makes sense to use the upper, almost horizontal, portion of the curve Kp, = f(Q). The range of this area may 
be extended either at the cost of increasing the coefficient A (whereby the gain coefficient is also increased) or 
of lowering ry (with simultaneous lowering of Kp: 


It should be mentioned that Conditions (6) and (27) for obtaining maximum power gain coefficients do not 
afford the simultaneous attainmentof minimum frequency distortion,since in this case the operative points on the 
curves Kp, = f(Q) would be only the point of inflection. 


SUMMARY 


One can not effectively determine the power gain coefficient of a magnetic amplifier destined for ac amp- 
lification by means of the apparent or operative power at the input. The existing formulas for this purpose do not 
allow one to judge the completeness of utilization of every amplifying device. Moreover, these formulas cannot 


be applied for every type of magnetic amplifier. 
It was suggested that the power gain coefficient be found as the ratio of the power developed in the load to 


the maximum power which the signal transducer matched to the load is capable of supplying. In an ideal magnetic 
amplifier controlled by some sinusoidal ac signal, maximum power gain (also maximum power in the load) is at- 


tained when the time constant of the amplifier control circuit equals the reciprocal of the signal's angular frequency, 


This condition determines the basic parameters of the amplifier, and allows a rational computation of them. If 

it is desired to obtain the maximum average power gain coefficient for a magnetic amplifier working in some fre- 
quency band from Q, to 9, then its time constant must be set equal to r = 1/9,Q,. If the amplifier's frequency 
distortion coefficient is given, its time constant can be determined from Equations (29) or (30). 
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CRITIQUE 


ON ESSENTIAL INADEQUACIES OF THE DEFINITIONS AND TERMINOLOGY OF AUTOMATIC CONTROL* 


A. L. Vertsaizer and V. M. lakovleyv 


(Daugavpils) 


Certain critical remarks are presented, bearing on the essential inadequacies 
of the definitions, terminology, and ambiguities presented by the structure of de- 
signs for systems of automatic control. 


The absence of unanimity onelements and principles, ambiguities and multiple meanings in the presenta - 
tion of designs, in the definitions and in the terminology of automatic control systems—all these are frequently to 
be found in contemporary literature, and particularly in textbooks, on the subject of automatic control. 


Let us turn to the textbooks on automatic control which have appeared during the last decade, The textbook 
authors named below adhere to a single system of definitions for automatic control systems, and for the nomen- 
clature of automatic controllers. However, the schematics of avtomatic control systems are given variously, not 
only in the sense of the sketches, but also, in our opinion, in the principles themselves. 


Consider book [1]. We read on page 19 that: 


Gout Qin "The controlled object, the regulating organ and 
fin the automatic controller are the elements of an 
, POb automatic control system. The automatic controller 
ol and the regulating organ together comprise the re- 


gulating device by means of which the controlling 
process is implemented. Certain authors consider 
the controlling organ as an element of the automatic 
L770 | controller. It appears to us that this is not com- 

12 pletely correct, since the controlling organ might 
Fig. 1. I is the automatic controller ba oct tate by 


Later, on page 32, we read: “A device which 
is used for the automatic displacement of a regulating organ in order to obtain a previously given nominal value 
of the controlled parameter is called an automatic controller." In Fig. 19 (page 44[1]) there is given a schematic 
of a controlled object with a linearly acting controller (Fig. 1), and in Fig. 29 (page 62[1]), there is given a 
schematic for a controlled object with a nonlinearly acting controller without feedback (Fig. 2). 


On the basis of the schematic of Fig. 1, M. V. Semenov draws the conclusion that ", . . a linearly acting 
automatic controller is comprised only of a sensing element." 


It should follow then from Fig. 2 that a nonlinearly acting controller consists of a sensing element and 
a servomechanism.. Thus, using M. V. Semenov's text, one may enunciate*the followings 


1. An automatic control system consists of three elements; regulated object, regulating organ and auto- 


* The editor considers very important the development of criticism by the readers in the direction indicated in 
this article. 


matic controller. The automatic controller and the regulating organ comprise the controlling device, i.e., besides 
the concept of automatic controller there is now introduced the concept of controlling device. 


2. The regulating organ is not an element of the automatic controller. 
M. A. Aizerman [2] says on page 23; "A system 
a of straight-line control consists of two elements, con- 
_ es Aone nected in a closed loop. One of the elements is the con- 


trolled object and the second is the controller. ... In 


| | the simplest cases, a nonlinear control system consists 
lg |, | of three elements: the controlled object, the sensing ele- 
iT | ment, reacting to deviations of the controlled quantity 

| a | from its required value, and a servomotor, controlled by 
| the sensing element.” In Fig. 7b [2] is shown the inter- 

| mi. | action schematic for a linear control system, and in Fig. 
| 


(renee | | a, 9b [2], the schematic for a non linear control system 
iy (Fig. 3 and 4). 


On page 22 [2] it is asserted that: “In a linear con- 
Fig. 2. Lis the automatic controller, trol system, the controller has one degree of freedom. 
II is the executive mechanism. It consists of one element which reacts to variations of 
the controlled quantity and immediately affects the re- 
gulating organ of the regulated object." As is clear from 
1 Fig. 3 and 4, M. A. Aizerman does not show in his sche- 
@ matic the action of the regulating organ as one of the 
elements of the automatic control system, as does M. V. 
Semenov. 


j In his book [3], E. P. Popov says on page 20; "A 


system of automatic control consists of a controlled ob- 
ject and a controller . . . By the regulator is meant the 

automatic device perceiving undesirable deviations of 

the controlled quantity x and so acting on the controlled 
; object that these deviations are removed... . For im- 

plementing its functions, the regulator first must have at 
its input a sensing element for perceiving the undesired variations of the controlled quantity and, second, the re- 
gulator must have at its output a controlling organ, exerting on the conttolled object a definite influence toward 
expunging the generated undesirable deviation of the controlled quantity. In addition to the sensing element (at 
the input) and the controlling organ (at the output), the controller may have various intermediate devices,” Cor- 
responding to this is given, in Fig. 8 (page 21[3]), a structural schematic (Fig. 5). 


Fig. 4. 1 is the servomotor, 2 is the sensing element, 
3 is the controlled object. 


In another book by E. P. Popov [4], Fig. 12 (page 23 [4]) shows the block schematic of an automatic control 
system with linear control (Fig. 6). Thus, according to E. P. Popov, a linear automatic controller consists of a 
sensing element and a controlling organ, but a nonlinear controller, in addition to the aforementioned elements, 
also includes intermediate devices. It is clear from just the three texts considered that there does not exist una- 
nimous agreement as to the elements comprising automatic control systems and automatic controllers, and that 
the block schematics given in these texts appear quite diverse. After reading the above-mentioned texts, one in- 
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voluntarily raises the question: to which portion of the controlling system should the controlling organ be referred? 
Is it an element of the automatic controller or of the controlled object, or is it necessary that it be considered an 
individual element of the system of automatic control? At firt glance, this question might appear trivial. How- 
ever, we feel that the given question is a cardinal one, for the following reasons: 


genre 
4 at 
Fig. 5. I is the controller; 1 is the Fig. 6. I is the automatic controller; 
disturbing influence f(t), 2 is the 1 is the disturbing influence, 2 the 
controlled object, 3 is the sensing controlled object, 3 the sensing ele- 
element, 4 is the controlling organ, ment, 4 the controlling organ, 5 the 
5 is the adjustment - y and 6 is the adjustment. 


set of intermediate devices. 


1) In studying automatic control systems, using the literature sources, one gets different opinions as to the 
basic elements of automatic control systems and, in some cases, as to the elements entering into the makeup of 
the automatic controller and the controlled object; 


Fig. 7. Fig. 8, 


2) The automatic regulators actually used in technology fall into two variants: one has built in both sensing 
element and controlling organ, and the other has only the sensing element. 


For example, the linear action, angular stress regulator, employed in aircraft technology, consists of a sens- 
ing element and a controlling organ, implemented in one box, while the regulator of the aircraft cabin temperature, 
@ nonlinear one, consists only, essentially,of a sensing element and intermediate devices. However, the first and 
the second both bear the name of automatic controllers (avtomaticheskogo reguliatora), 


In our opinion, it is necessary to consider as automatic controllers devices consisting of the following elements: 
a sensing element, intermediate devices and a controlling organ (reguliruiushchii organ). In support of this asser- 
tion, we consider a system for maintaining rotational stability of a heat engine. 


With manual regulation (Fig. 7) a human observes the rotation of the engine by means of the indicator P, 
compares the actual rotation with the desired value, and by means of a valve Z increases or decreases the supply 
of energy to the engine, D. This is not automatic regulation, but the human, together with the enumerated devices, 
does execute the role of automatic controller. To make the system automatic, it is necessary to set up a centrifu- 
gal mechanism and connect this to the valve. 


| 
| 


A human alone does not play the role of an automatic controller. He needs supplementary elements: a 
counter of revolutions and a valve. In automatic regulation by the centrifugal mechanism CM (Fig. 8), one part 
of the automatic controller does not appear. For this, it is necessary to have a connection to che valve via a ki- 
nematic linkage, L. 


Consequently, it is correct to state that a linear automatic controller consists of two elements, a sensing ele - 
ment and a controliing organ, but a nonlinear automatic controller consists of a sensing element, intermediate 
devices and a controlling organ. In.our opinion, attention should also be paid to terminology. In scanning block 
schematics of automatic control systems, we convinced ourselves that one and the same element of an automatic 
control system, for example, the sensing element (chuvstvitel'nyi element), is denoted differently in different 
books. E. P. Popov and M. V. Semenov use the term “sensing element," V. A. Bodner in [5] in Fig. 2.2 (page 29 
[5]}) calls this same element “measuring sensing element of the controller", and in his second edition [6], in Fig. 
11.3 (page 31 [6]), calls this same element “measuring device”. 


A. A. Voronov in [7] calls this same element the “measuring organ". One and the same element of the system 
is called “controlling organ" by A. A. Voronoy [7], “controlling element" by E. P. Povov [3], M. V. Semenov [1] 
and K. V. Egorov [8], and "amplifier" by V. A. Bodner[5]. The driving device is called "seryomotor" by V. A. 
Bodner and M. V. Semenov, “power element” by K. V. Egorov, and “drive” by E. P. Popov. It seems to us that 
enough has beeen said in this connection to warrant drawing certain conclusions. It is time to develop unified 
concepts and definitions, and to adhere to one terminology in questions of automatic control. The question is 
placed open to general discussion. 


It would seem that if the authors of the texts mentioned were to enter into this discussion, and whole -heart- 

edly, the problem would be solved. 
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CHRONIC LE 


CONFERENCE ON ELECTRICAL CONTACTS 


During November 26-28, 1956, a conference on electrical contacts was held in Moscow. The conference 
was organized by the IAT (Institute for Automation and Remote Control) of the AN SSSR, in conjunction with 
various research institutions. One hundred seventy people from 53 organizations participated: 33 research institues 
and 20 plants turning out electrical apparatus and contacts. At the conference 25 papers were heard, after which 
a 22-man panel discussion was held. At its conclusion, the conference adopted several resolutions, 


The aims of the conference were to clarify the present state of research work in the domain of contacts and 
the phenomena related to contact functioning, the investigations into increasing the useful life of contacts and de- 
termining the reliability of their use, and the questions of production and utilization of contact materials. 


Electrical contacts figure among the most wide-spread and critical of the elements in electrical designs, On 
contact properties depend, in large measure, the reliability and length of useful life of electrical apparatus. It is 
in this connection that, recently, there has appeared a large quantity of work in the determination of useful life 
and domain of rational employment of contact materials and various contact types, as well as work investigating 


the phenomena connected with electrical contacting functioning, in so far as those processes were hitherto not 
completely clear. ' 


In the conference, the first of its kind in the Soviet Union, there was considered a wide range of questions, 
having to do both with electrical contacts directly and with related areas, 


The conference papers were devoted to the following questions. 
1. The physics of the processes entering into the contact function. 
2. Contact computations, and methods for testing contact apparatus. 


3. Current contact materials, their characteristics, properties, production, and recommendations for their 
usage under various conditions. 


In the papers devoted to the physics of the processes entering into contact functioning, there were explained 
several endeavors in the investigation of spark discharge and the electric spark machining of metals. 


S. L. Mandel'shtam, V. P. Shabanskii and N. K.Sukhodrey (Moscow) showed that spark discharges are charac- 


terized by high energy density produced at the electrodes, whereupon the leakage from the cathode occurs explo- 
sively. 


L. S. Palatnik (Kharkov State University) showed that the fundamental role in the transfer phenomenon in 
electric spark machining is played by the physico-chemical nature of the electrodes and their environment, and | 
stated criteria for the transfer and the mutual interaction of the electrodes. f 


B. N. Zolotykh (Central Research Laboratory for Electric Spark Machining of Materials, AN SSSR) showed 
that the amount of erosion is related to the heat constant of the material, and proportional to the energy supplied 
in the pulses. 


N. V. Afanas'ev (White Russian Polytechnic Institute) compared the physical constants of materials with their 
erosion stability, and attempted to establish certain connections between the magnitude of the erosion and the heat 
of evaporation. In this group should also be includedgthe papers of I, V. Kragel'skii (Institut Mashinovedeniia AN 
SSSR), who demonstrated both theoretically and on models the actual dependence of the area of the contact surface 


| 


on the properties of the surfaces the paper of A. I. Pugin (Metallurgical Institute AN SSSR), investigating the character 
of the initial formation of the contact, the potential difference on the contact, and the temperature in the neigh- 
borhood of the contact during a welding process; and the paper of M. A. Razumikhin (Moscow) which presented 
methods of coping with bridge and spark erosion, and on the influence of the environmental pressure on the func- 


tioning of a contact. 


Among the papers dealing with the second group of questions should be included the paper of B. S. Sotskov 
(IAT AN SSSR) which expounded an approximate method for computing relay contacts for different conditions of 
use, taking into account their erosion and the probability of their breakdown as a function of their length of ser- 
vice, the paper of O. B. Bron (Leningrad), which gave an analysis of the conditions of functioning of contacts in 
contactors and in automatic cut-outs, the dimensions and parameters of commutating links with and without arc 
quenching, and the papers of G. O. Feiler (Moscow) and V. I, Fiks and M, A. Meerovich (Plant ATE-1 Ministry of 
the Automotive Industry), presenting the results of industrial tests of different contact materials, and the paper of 
A. V. Gordon (Moscow), reporting on the requirements presented by contacts for aircraft contactors, and stating 
methods for accelerating contact testing. 


The third group of questions was illuminated by a large number of papers dealing with new contact materials, 
with certain peculiarities of previously known materials, and with the replacement by these of the precious metals. 


A. B. Altman, I. P, Melashenko and E. S, Bustrova (Moscow) reported on existing methods of producing me- 
tallo-ceramic contacts, and showed that during their operation the surface layer of metallo-ceramic contacts 
was changed in structure and in phase composition. 


I, N. Frantsevich and O. K. Teodorovich (Metallurgical Institute AN USSR) state that the most economical 
and general method of producing metallo-ceramic contacts is the impregnation method, and that for obtaining 
material with the very best characteristics it is necessary to compress the composition material. 


I. E, Dekabrun (IAT AN SSSR) explained the results of work in the determination of certain characteristics 
of metallo-ceramic contacts, and recommended certain usages of such contacts in the domain of weak currents, 


V. V. Usov and E. M. Murav"eva (Moscow) reported on contact materials obtained by internal oxidation of 
certain alloys, the properties of which approximate to the properties of metallo-ceramic contacts, but the tech- 
nology of preparing which are significantly simpler and cheaper. 


A. B. Al‘tman and E. S. Bustrova (Moscow) reported on the possibility of using aging metallurgical alloys 
for contacts, and showed that their stability under wear is very high. o 


V. Y. Usov and M. D. Povolotskaia (Moscow) presented data on the reasons for corrosion of tungsten contacts; 
E. A. Shumskaia(NII Avtopriborov Ministerstva avtomobil ‘noi promyshlennosti) showed that the erosion of tungsten 
contacts depends on the deformation of their atomic lattices. 


A. A. Rudnitskii (Metallurgical Institute AN SSSR) presented data on the thermoelectromotive force of a 
series of precious metals and their alloys, which plays a large part in the behavior of low current contacts. 


E. S. Kirillova (Moscow) compared the characteristics of palladium-silver and other alloys, determined with 
the object of using these alloys for sliding contacts with very small pressures. 


E. K. luferova (Moscow) made recommendations on the use of a series of new materials for sliding contacts 
in synchronous transmission systems. 


V. A. Mitiushey and V. V. Kozlov (Moscow) made a survey of the research work being done in industrial 
plants, and reported on the state of contact and contact material production in these plants. 


The discussions brought out the great value of the conference for the wide interchange of experience and the 
dissemination of results obtained. Moreovet, a series of remarks was made in the papers, reporting on a host of new 
labors, carried out in the area of contact research, and a number of assertions were aired. 


Among the resolutions adopted by the conference should be mentioned the necessity for periodic sessions (one 
every two or three years) of the conference on electric contacts, and also the necessity of coordinating research 
work in the investigation of processes occuring in contacts, in the investigation of useful life, in the development 
of new contact materials, and in the establishement of standard methods for testing contacts. 


I, E. Dekabrun 


